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PREFACE 

This  reference  document  was  written  by  Mr.  Gilbert  C. 
Binninger,  Dr.  Paul  J.  Castleberry.  Jr.  and  Ms  Patsy  M.  McGrady 
under  the  supervision  of  Mr.  John  W.  Burfening  of  tne  Physical 
Vulnerability  (Nuclear)  Branch  of  the  Targets  Division,  Defense 
Intelligence  Agency.  This  document  presents  in  detail  the 
mathematics  of  selected  portions  of  the  Defense  Intelligence 
Agency  publication  AP-550-1-2-69-INT,  "Physical  Vulnerability 
Handbook  - Nuclear  Weapons  (U),"  1 June  1969.  The  appendices 
present  the  logic  flow,  necessary  equations  and  constants 
needed  for  computerization  of  the  PV  Handbook.  Illustrative 
programs  written  in  BASIC  language  are  included  with  each 
appendix. 

Suggestions  for  additions  or  improvements  to  this  document 
should  be  transmitted  to 
Director 

Defense  Intelligence  Agency 
ATTN:  DI-7D 

Washington,  D.C.  20301 
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The  Defense  Intelligence  Agency  publishes  "Physical 

Vulnerability  Handbook  - Nuclear  Weapons  (U),"  AP-550-1-2-INT, 

# 

for  the  use  of  operational  planners,  target  officers,  and 
physical  vulnerability  analysts  concerned  with  the  employment 
of  nuclear  weapons  and  the  prediction  of  their  effects.  The 
methodology  presented  in  the  PV  Handbook  allows  for  target 
hardness,  weapon  yield,  weapon  delivery  system  accuracy, 
height-of -burst,  and  doctrinal  requirements  to  be  considered 
in  practical  targeting  problems. 

The  ancestors  of  the  PV  Handbook1 * 2 omitted  mathematical 
formulations  to  decrease  their  bulk  without  detracting  from 
their  day-to-day  effectiveness.  The  result  is  that  today  the 
"Physical  Vulnerability  Handbook  - Nuclear  Weapons  (U) " presents 
only  tabulated  data,  graphs , and  other  computational  aids  for 
solution  of  practical  problems  and  almost  completely  omits 
the  mathematical  formulations  needed  to  understand  the  system. 


"Target  Analysis  for  Atomic  Weapons",  AF-62B202,  by  the 
Physical  Vulnerability  Division,  Directorate  for  Intelligence, 
Headquarters  United  States  Air  Force,  PV  TM-14,  30  June  1954. 

2"Nuclear  Weapons  Employment  Handbook  (U),"  AFM  200-8,  by  the 
Physical  Vulnerability  Division,  Director  for  Targets,  Assistant 
Chief  of  Staff,  Intelligence,  Headquarters  United  States  Air 
Force,  1 May  1958,  SECRET/FRD.  ‘ 


1 


This  publication  is  intended  to  be  a reference  work  con- 
taining the  detailed  mathematical  formulations  whifh.  constitute 
the  groundwork  of  the  Physical  Vulnerability  system  and  in 
particular  to  assist  those  doing  computer  calculations.  This 
presentation  is  especially  needed  to  document . Change  1 to  the 
PV  Handbook  which  was  issued  1 September  1972.  Change  1 
dropped  the  circular  coverage  function  and  adopted  the 
cumulative  lognormal  function  to  describe  probability  of 
damage  as  a function  of  distance  from  the  detonation  point. 
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II.  ANALYTIC  APPROXIMATION  TO  THE  DISTANCE  DAMAGE  FUNCTION 
A.  Factors  Leading  to  Choice 

In  1951  the  Hiroshima  and  Nagasaki  data  were  analyzed3 
to  determine  the  damaging  overpressures  and  variation  in  damage 
with  respect  to  overpressure  for  the  various  types  of  structures 
found  at  Hiroshima  and  Nagasaki. 

These  data  were  insufficient  to  determine  the  precise 
mathematical  form  of  the  probability  of  damage  as  a function  of 
pressure.  However,  after  the  evaluation  of  several  analytic 
approximations , the  cumulative  lognormal  function  (a  cumulative 
normal  function  with  logarithmic  variable)  was  chosen  as  the 
best  fit  to  the  collected  probability  of  damage  versus  pressure 
data.  * 

Probability  of  damage  versus  range  curves,  referred 
to  as  distance  damage  functions,  are  generated  for  any  desired 
height-of -burst  (HOB)  by  combining  the  pressure-damage  curves 
with  pressure  vs.  range  curves  for  the  desired  heights-of-burst 
(HOB) . These  distance  damage  functions  are  numerically  inte- 
grated to  obtain  the  weapon  radii  for  Part  I of  AP-550-1-2-69- 
INT  (See  Appendix  A) . 


3"A  Classification  of  Structures  Based  on  Vulnerability  to 
Blast  from  Atomic  Bombs",  Technical  Memorandum  No.  4,  Physical 
Vulnerability  Eranch,  Air  Targets  Division,  Directorate  for 
Intelligence,  United  States  Air  Force,  2 March  1951. 

3.  e#  * ‘ • . < » « • 
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To  simplify  the  probability  of  damage  calculations 
in  Part  IV  of  AP-550-1-2-69-INT  when  weapon  delivery  error  is 
included,  the  actual  distance  damage  function  described  above 
is  approximated  by  an  analytic  function.  In  order  to  give 
sufficiently  accurate  results  for  practical  applications, 
this  analytic  function  must  comply  with  certain  characteristics 
of  the  actual  distance  damage  function.  PV  TM-144  showed  that 
the  chosen  mathematical  model  must  closely  approximate  the 
actual  distance  damage  function  in  the  maximum  level  of  damage 
expected  in  the  neighborhood  of  the  ground  zero  and  in  the 
distance  at  which  there  is  a 50%  probability  of  damage.  Also, 
the  relative  slope  and  general  shape  of  the  analytic  function 
must  agree  within  reasonable  limits  with  the  slope  and  shape 
cf  the  a stool  distance  damage  c ..rvc . 

Historically  the  Circular  Coverage  Function  (CCF)  has 
been  used  for  this  approximation.  However,  for  targets  havener 
damage  sigmas  greater  than  0.3C  the  CCF  :.c  ur-obt-  to  satisfy 
the  requirement  that  probability  of  damage  be  near  100%  at 
ground  zero.  Therefore,  beginning  with  Change  1 to  the  1969 
PV  Handbook,  the  complement  of  the  cumulative  lognormal  function 
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is  used  as  the  analytic  approximation  to  the  distance  damage 
function.  These  two  approximations  for  the  distance  damage 
function  are  virtually  identical  for  %-  JO  , differ  at  most 
4 percentage  points  for  0j  = .20  , and  differ  at  most  6 percentage 
points  for  0j».3O  . These  differences  do  not  cause  significant 
differences  in  calculated  probabilities  of  damage  when  weapon 
delivery  error  is  included.  However,  significant  differences 
in  calculated  probabilities  of  damage  can  occur  for  Cj’s  of  .40 
and  .50. 


B.  Analytic  Form  of  the  Lognormal  Distance  Damage  Function 
In  Part  IV  of  AP-550-1-2-69-INT  the  actual  distance 
damage  function  is  approximated  by  the  complement  of  the 
cumulative  lognormal  function  in  order  to  simplify  probability 
of  damage  calculations.  The  lognormal  density  function  with 
variable  r and  parameters  cL  and  p is  expressed  mathematically 
as 


for  r > 0 where  At  is  the  natural  logarithm.  The  standard 
notation  has  r and  CL  interchanged  in  the  exponent,  however, 
the  forms  are  equivalent.  This  density  function  is  plotted  in 
Figure  1 below  for  several  values  of  3 for  c £.  = / . 


FIGURE  1, 


over 


It  can  be  shown  that  cL(^  ) is  a density  function 

the  domain  (0,oo)  by  letting  - -Jp  Jlr\  ( ^ ) and  noting  that 


ridi  **(£>]=  -«  and  rli”  [t u ^)] = 


cO 


Therefore , 


j-p(rj0.,?)dr  ~j^±=,  e dtf  = jygf  e ' 


The  parameter  o£  , which  is  the  median  of  the  lognormal 

► 

density  function,  is  r , the  distance  from  ground  zero  at 

which  there  is  a 50%  chance  of  achieving  a specified  level  of 

A. 

damage,  i.e.,  j oc  p)dr  = Vz  • The  parameter  ^3  is  the 
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t 

I 


s* 


standard  deviation  of  Jyi r . 

The  cumulative  lognormal  function  (i.e.,  the 
distribution  function)  is  expressed  as 

Pfr)  =-  I d?*j-6irC  ^ 

° Ur) 

where  H (0  = ■Jj“  /tl  (” r-)  ' 

A graphical  representation  of  P (r)  is  shown  in  figure  2. 


FIGURE  2 


Since  the  distance  damage  function  Pd(r)  is  monotonically 
decreasing  it  is  expressed  as  the  complement  of  P(r)r. 


Pd(r)  = l-P(r) 


P,(r)  = | 


e 

Z(r) 


-n% 


* t(r) 

dt  = j"vfr  e 


•t 


-OQ 


where  z(r) 


The  distance  damage  function,  P/r)  , is  depicted  in  figure  3. 

I.o 

P/O 


r — * 

FIGURE  3 


Note  that  the  damage  function  P/r)  is  the  complement  of  a 
cumulative  distribution  function,  not  a density  function,  i.e., 


o 


Since  Pd  (r)  is  a cumulative  function,  it  represents  the 
probability  that  a target  will  receive  at  least  a specified 
level  of  damage.  For  instance,  if  P/r)  is  being  used  to 
describe  severe  damage,  1-P/r)  is  not  the  probability  of 
survival,  but  the  probability  of  achieving  less  than  severe 
damage . 
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III.  PARAMETERS  USED  TO  DESCRIBE  THE  DISTANCE  DAMAGE  FUNCTION 

A.  Definitions  of  WR  and  , Parameters  of  the  Distance 
Damage  Function 

Given  the  analytic  form  of  Pd(r)  from  Section  II. B., 

the  parameters  oL  and  |3  uniquely  specify  the  damage  function. 

Historically,  however,  the  DIA  PV  methodology  has  used  the  two 

quantities  weapon  radius  (WR)  and  distance  damage  sigma  ( ) 

insteaJ  of  oL  and  fi  to  describe  Pd(r).  The  quantities  WR  and 

are  constructed  from  the  first  two  moments  about  the  origin 

of  the  lognormal  density  function.  Specifically,  WR=  V< »“*>’ 

and  cn2  * — . , r where  the  first  moment  <r>  and  second 

< r s 

moment  <r*>  are  defined  below. 


where  the  first  moment  <r>  and  second 


<r>  s = JP/r)  dr 

;r*>  ■ j2rP*(r)d' 


The  last  transition  on  each  line  is  achieved  through  integration 
by  parts.  The  procedure  for  the  numerical  calculation  of  <r> 
and  <r4>  for  any  distance  damage  function  is  presented  in  appendix 


The  functional  relationship  between  oL  and  |3  and  WR 


•w'.yr-vmni.  ,n."»u,i 


m»y,u 


and  Oj  will  be  derived  in  part  C below. 

B.  Conceptual  Descriptions  of  WR  and  <7^ 

The  weapon  radius  is  the  square  root  of  < r*>  , the 
second  moment  of  the  density  function  associated  with  the 
distance  damage  function.  A more  conceptual  definition  states: 
Given  a uniform  distribution  of  like  targets,  the  WR  is  the 
radius  of  a circle  centered  at  the  GZ  that  contains  as  many 
targets  undamaged  to  a specified  level  inside  as  there  are 
targets  damaged  to  a specified  level  outside.  Thus  if  the 
undamaged  targets  inside  the  circle  are  replaced  with  the 
damaged  targets  outside  the  circle,  the  circle  of  radius  WR 
would  enclose  an  area  entirely  damaged  to  the  specified  level. 
Expressed  mathematically, 

ESX  ™ama9ed  - J / (I  - P/r>)  r dr  36 

o o 

2V  « 

Number  damaged  _ ( f Q / \ . iA 

outside  WR  ' J J rd(r'  rdr  dy 

0 WR. 

Equating  these  and  simplifying, 

to  oo 

WR2  = j Z r fjco  clr  = | rZ  f.(r-  cl}  p)  dr  = <ra> 

o o 

i.e.,  WR2  is  equal  to  the  second  moment  about  the  origin  of  the 
lognormal  density  function. 

The  mean  area  of  effectiveness  (M.A.E.)  of  a weapon 


* 


10 


is  defined  as  the  area  over  which  a weapon  on  the  average 
achieves  at  least  a specified  level  of  damage: 
iir  to 

m.a.e.  “ £ f PjO)  rdrde  =*  IT  WR2 

o o 

The  mean  area  of  effectiveness  of  a nuclear  weapon  is  a circle 
of  radius  WR. 

The  second  basic  parameter  of  the  PV  system  is  the 
distance  damage  sigma,  0^  , a number  between  0 and  1. 

<ra> - <r>z  _ <r>* 

°d  " <r*>  1 “ <r*> 

CT^Z  is  a measure  of  the  variance  of  the  density  function.  It 

is  the  variance  made  dimensionless  by  division  with  the  second 

moment.  Small  o^* s indicate  a relatively  rapid  fall  off  of  the 

damage  function.  Large  0^*5  indicate  a more  gradual  fall  off. 

The  damage  sigma  should  never  be  confused  with  the  derivative 

(slope)  of  the  damage  function. 

C.  Evaluation  of  WR  and  <3^  for  the  Lognormal  Distance 
Damage  Function 

In  order  to  evaluate  WR  and  in  terns  of  the  para- 
meters CX-  and  of  the  lognormal  distance  damage  function, 
it  is  necessary  to  evaluate  the  first  two  moments  of  the 
lognormal  density  function.  ■ . ■' 

11 

»> 


The  second  moment, <r*>,  of  the  lognormal  density  functior 


is  defined  as 


«D 

<r*>  a j ' rz  -per,  CL,  p)  dr 


Again  using  the  substitution  r 


~SO  “ t 


this  becomes 


<rx> 


-f 


VztP 


. dlJ 


By  completing  the  square. 


i r, 

W 


« 

Ox>  = rso  e2**  ( 


YztT  e 


Since 


1 -x(y+*p)' 


Ju  "=*  l • the  second  moment  is 


< r*>  = r*  e*P 


Using  the  above  results  f or  <r*>  and  the  definition 
WR  = Y < rl>'  , the  WR  is  given  by 


I3  I 
= r,a  e1  * Ot  e' 


" 1 ^ ^"rvr-pzri^* g^n-»  ■ztmrrrr&'rr  ‘-1-  ».v« 


The  above  results  for  <r>  and  <c7->  can  be  Substituted 
into  the  equation  for  (J^z  to  obtain 


<r*>  -<r>x  ._ 
<rl> 


This  can  be  inverted  to  obtain  |5  for  the  lognormal  distance 
damage  function  if  the  distance  damage  sigma  is  known.  The 
inversion  gives 


yS  = V-"/n  ( l-<rf) 

The  expressions  for  WR  and  0^Z  can  be  combined  to  give  r$Q  = ot 
if  WR  and  cq  are  known.  Specifically," 


rso  = cc  = wr  ( i - O 
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IV.  COMBIK-D  weapon  effects 

When  a target  is  primarily  susceptible  to  one  weapon 
effect,  then  the  WR  and  OJ  of  that  effect  alone  are  used  in 
predicting  damage  to  the  target.  For  example,  a single  story 
wood  frame  dwelling  is  primarily  susceptible  to  overpressure 
and  therefore  it  is  designated  a P (for  overpressure)  target. 

In  this  case  the  probability  of  damage  to  the  target  is  cal- 
culated using  only  the  appropriate  WR  and  oj  for  overpressure. 

In  some  cases,  however,  two  or  more  weapon  effects  may 
significantly  contribute  toward  damaging  the  target.  For 
example,  personnel  casualties  may  result  from  blast  effects 
and/or  radiation  effects.  In  these  situations  the  weapon 
radius  has  a larger  value  than  it  would  have  for  either  effect 
alone.  The  DIA  methodology  assumes  the  independence  of  effects, 
Therefore,  the  combined  effects  weapon  radius  is  estimated 
without  allowance  for  the  possible  greater  susceptibility  of 
a target  to  one  effect  because  of  exposure  to  another. 

If  three  independent  effects  (i  = 1,  2,  3)  are  considered, 
and  the  damage  function  for  each  is  given  by 


„ rir]  -1*4; 
* ) Tin  e d1 


where 


«« • * *.  CM1) 
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then  the  combined  distance  damage  function  is 


Vr>  - l-0-PJ,(r))(l-R1,(r)Xl-PJ,(r)) 

- I * f I -[  l -(  I - Pj(r)X  I"  Pj  «)]][  I - Pj/D) 

- l-(l-f51>lfr)Xl-P41<p)) 

wherepd12(r)  = %<rd . 

The  above  formulation  represents  the  combining  of  two  of  the 
effects  and  then  combining  that  resultant  effect  with  the 
third. 

The  weapon  radius  of  two  combined  effects,  WRJ2,  is5 

co  m 

WR*Z  = j"  2 r Pd#i( r)  dr  = WRf  + wR*  - 2 | r /j  <r)  P^CO  d r 

O O 

wr*  = wi?x  Bfz'fw^e*))  + ml  B(2(mxePx)')  (equ.  d 

i 

i 

where  e • c 

r i 

b(h»)  = j me  Ji 

<(r) 

5The  derivations  cf  cguations  1 ar.d  2 are  presented  in  detail 
in  Appendix  g. 
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and 


/ , ' i , ( wrl  e’*1 


The  first  moment  when  two  independent  events  are  combined, 
<ria>  , is  derived  in  the  same  manner. 


•* 


<f'ia>  - <f,>B«(WR,))  + <rz>  B(2,'(wrJ) 


(EOU . 2) 


The  distance  damage  sigma  for  two  combined  effects,  Oj  , is 

1 3. 

found  using 


A computer  program  for  calculating  the  weapon  radius 
and  damage  sigma  for  two  combined  effects  is  presented  in 
Appendix  C . c V. 
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V.  WEAPON  DELIVERY  ERROR  ■ . v % ^ = | 

The  circular  error  probable  (CEP)  is  a measure  of  weapon 
system  accuracy.  It  is  the  radius  of  a circle  centered  at  the 
desired  ground  zero  (DGZ)  within  which  50%  of  the  impact  points 
will  fall  if  the  distribution  of  impact  points  is  assumed  to  be 
normally  distributed  about  the  DGZ.  Expressed  mathematically: 

JOT  CEP  _r}/  z 

x=  ( ( e rJrde 


The  standard  deviation.  O',  of  the  impact  point  distri- 
bution is  often  expressed  in  terms  of  the  CEP  as  follows: 


VI 


PROBABILITY  OF  DAMAGE  CALCULATIONS 


A.  Point  Targets 

By  using  the  damage  function  presented  in  Section  II.  B, 
it  is  possible  to  predict  the  probability  of  damage  to  a point 
target  if  the  offset  distance  from  the  detonation  point  to  the 
target  is  known.  In  general,  however,  all  that  is  known  is  the 
distance  from  the  aim  point  to  the  target.  The  actual  distance 
from  the  weapon  detonation  point  to  the  target  is  uncertain  due 
to  inaccuracies  in  the  weapon  delivery  system.  Thus  the  results 
of  an  attack  can  not  be  predicted  with  certainty;  all  that  can 
be  predicted  in  most  cases  is  what  is  most  likely  to  happen,  or 
what  will  happen  on  the  average.  The  PV  system  calculates  an 
average  probability  of  damage  by  weighting  the  probability  of 
damage  for  each  possible  detonation  point  by  the  probability 
that  the  weapon  lands  at  that  detonation  point.  The  probability 
curves  in  Part  IV  of  AP-550-1-2-69-INT  were  derived  assuming 
that  the  distribution  of  actual  detonation  points  about  the  aim 
point  is  described  by  the  circular  normal  distribution,  i.e., 
the  probability  that  a weapon  lands  in  a small  area,  £A,  a 
distance  O from  the  DGZ  is  given  by 


I 

2rr<rx 


e 


A A 


where  a* is  the  standard  deviation  of  the  circular  normal 


o 


point  target  located  a distance  x from  the  DGZ, 


P*  Z P»  • 


•n 

*r««. 


I 

2TTcra 


-rx 


20" 


AA 


2ir  » 


P* 


(r*+-x*-2rx  cose) 


o a 


To  calculate  P in  Part  IV  of  AP-550-1-2-69-INT,  Pd(r) 
is  approximated  by  the  complement  of  the  cumulative  lognormal 
function  presented  in  Section  II. B. 


*r)  u*/ 

C 1 

- \ -we  ^ 


where  z(r)  - -jj-  /n  (•£*)  * ^ . 

The  value  of  Pd(r)  is  calculated  at  increments  of  r 
by  use  of  the  error  function  (erf)  definition 


Pd(r)  "■  *5  + *5  HH  er^(^) 


where  z = “ 

The  erf  can  be  approximated  by  erf  (u)  = '-(1) 

k 

where  u > 0 , D = 1 •+»  6^  U. 

i«l  4 < ‘ ‘ ca 
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and 


0.0705230784 


ei  = 

e2  = 0.0422820123 
e j = 0.0092705272 
e4  = 0.0001520143 
e s = 0.0002765672 
e6  = 0.0000430638 


I .Ztt  (V*+  X*-2r*  cose) 
The  term  x j e 


rdrd©  is  found  by 
using  series  expansions6  of  the  zeroth  order  modified  Bessel 


I /rirr  rx 

function  of  the  first  kind,  exp(=-~  Cos  0)  d6 , then 


multiplying  by  r exp  (- — - (r2  + x2)). 

20 

,2ir 


rx 

If  - 3.75  then 


i r 


-~\  exp  (~  cos  0)  d9  = 1 + 3.5156229*jt  + 3 .0899424* 2 

' + 1.2067492*j J3  + 0.2659732*jj4 

+ 0.0360768*j1s  + 0.0045813*j j6  +£ 

where  £.  < 1.6  x 10" 7 and  j,  -( — — — A 

\ 3 .75a2y  * 


_ - rx  • * • * , 

if  > 3.75  then  & ' ' 1 -V 

^ $2\xP(&coSe)cie=(0)  exp(p)  [0.39394228  + 

' . ■ + 0.01328592*j2 + 0.00225319*j 2 2 

- 0 . 00157565*j  2 3 + 0 . 00916281* j 2 4 

- 0 . 02057706* j 2 5 + 0 . 02635537* j 2 ( 

- 0 .01647633* j 2 7 + 0 .00392377* j 2 8 + 


6E.  E.  Allen,  Analytical  Approximation,  Math,  Tables  Aids  Comp 

8,  240-241  (1954) . 
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where  £ < 1.9  x 10”7  and  j2  = 


2.1  So1 
r*x 


The  original  integral  for  the  average  probability  of 
damage,  P,  may  be  thought  of  in  the  following  manner. 


Zrr  oo 


2TTCT 


->  ( ( PdM 


r xz-  zm  cos& 
2<r* 


D 

rJrd9- 


where  a and  b are  selected  such  that  when  r < a,  f (r)  s 0 and 
when  r > b,  f(r)  z 0 . The  integral^(r)dr  can  be  evaluated 
using  Gauss-Legendre  quadrature.  Using  this  method: 


( , . . (b-cO  WYb-a)  2l -tk +-<M 

J t(r)Ar  = i L 2*  ) £ 


where  w^  and  have  been  derived  for  various  n and  £ is  the 
accumulated  error.  Use  of  n=10  will  result  in  £ £ .0001. 

Appendix  E contains  a computer  program  based  on  the 
above  formulation  for  calculating  the  probability  of  damage  to 
a point  target  at  any  distance  from  the  DGZ.  c 

Appendix  F presents  a computer  program  for  calculating 
the  required  offset  distance  from  the  DGZ  to  the  target  to 
achieve  a specified  probability  of  damage. 


V'**' yt iv  '&-1--1 " m.1  Ij  ijiibw w^jj^wyip^^ mt# 
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B.  Area  Targets  . t „ ° xr 

Average  p.  cbabilities  of  damage  to  area  targets  are 

obtained  by  dividing  the  area  target  into  small  cells  which 
can  be  treated  as  point  targets,  calculating  the  probability  of 
damage  to  each  cell,  weighting  these  probabilities  by  the  area 
of  the  cell  or  the  portion  of  the  target  in  the  cell,  and 
averaging  the  results . This  approach  is  discussed  in  the 
following  paragraphs  for  normally  distributed  and  uniformly 
distributed  area  targets. 

1.  Circular  Normal  Distribution 

Some  area  targets,  such  as  population  centers, 

exhibit  a concentration  of  target  elements  in  the  center  which 

'■  '«  < > 

tends  to  become  less  as  the  distance  from  the  target  center 
increases.  The  distribution  of  target  elements  in  this  case 
can  be  well  described  by  a circular  normal  distribution.  A 
P-95  circle  is  used  in  Part  IV  of  AP-550-1-2-69-INT  to  describe 
a normally  distributed  population.  The  P-95  is  the  radius  of 
the  smallest  circle  which  encompasses  at  least  95  percent  of 
the  population  being  considered.  If  the  target  is  not  a pop- 
ulation target,  the  95%  radius  may  be  referred  to  as  an  R-95 
or  other  equivalent  term.  Expressed  mathematically 


.15 


ZTT  P-1S  _ _tl 

f f d?  tjtde 


This  method  of  calculating  the  probability  is  quite 
slow,  however.  A faster,  preferred  method  involves  combining 
the  target  and  weapon  delivery  distributions  into  a joint 
distribution  and  then  calculating  the  probability  as  if  for  a 
point  target  using  the  CEP  of  the  joint  distribution.  This 
approach  is  possible  because  both  the  delivery  error,  f(r), 
and  the  target  density,  f(t),  are  normally  distributed  and 
thus  their  joint  density  function  f(r,t)  is  also  normally 
distributed.  The  variance  of  a joint  distribution  of  inde- 
pendent random  variables  is  the  sum  of  the  variances  of  those 
random  variables.  Since  the  delivery  error  and  target  density 
are  independent,  f(r,t)  = f(r)  * f(t),  then  b£“tr+0^*'  where 

X ' . t cf  x 

<7^  describes  the  variance  of  the  joint  distribution,  <T  des- 

c r r-  X " 

cribes  the  variance  of  the  delivery  error,  and  0^  describes 
the  variance  of  the  target  density.  Therefore,  the  CEP  of  the 


'wi-n’i  Mnywff‘1— ii'Wi;1  y»'P  v '* 


w«yg  *wpf»gp  ■^■y^iupwf  >y ■■  r'n* 
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joint  distribution,  the  adjusted  CEP  (CEPa) , is  given  by 


■%r 


CEP  = 1. 1774  cr  = V Cepx  t-  .231  (P-^s)1 

CL  CL 

where  .ZSI  J = (|.|774  o;)*  is  the 

target  CEP2. 

A computer  program  based  on  this  formulation  for  cal- 
culating probabilities  of  damage  to  normally  distributed  area 
targets  is  presented  in  Appendix  E. 

2.  Uniform  Distribution  ' 

For  uniformly  distributed  area  targets,  the  method 
used  to  calculate  precise  probabilities  of  damage  again  begins 
with  the  division  of  the  area  target  into  small  cells.  All  of 
j the  cells  that  are  equidistant  from  the  DGZ  will  have  equal 

\ probabilities  of  damage.  Referring  to  figure  5,  if  concentric 

circles  are  drawn  about  the  DGZ  with  radius  r^  where  ri  = ri_^  + 
Ar  and  where  Ar  is  small,  then  the  target  cells  located  in  the 


If  is  the  area  of  intersection  of  the  i—  annulus  and  the 
target,  then  the  average  probability  of  damage  to  the  uniform 


target  is 


2 p.«; 

where  the  sum  of  the  a^  is  equal  to  ir(TR)2  and  TR  is  the  tar- 
get radius. 

A computer  program  for  calculating  probabilities  of 
damage  to  uniformly  distributed  area  targets  based  on  the 
above  method  is  presented  in  Appendix  G. 

Good  approximations  to  probabilities  of  damage  to 
uniformly  distributed  area  targets  can  be  obtained  by  treating 
these  targets  as  point  targets  and  replacing  the  weapon  delivery 
CEP  with  an  adjusted  CEP,  CEPa,  where 

CEPa  = 1CEPZ+  A TR*  for  + X +CEP) 


CEPo.  * ~JCEP*+  Tr1  for  TP>(uJfZ  + * -hCEP)  . 


These  CEP  adjustments  were  empirically  derived.  They  have  no 
basis  in  mathematical  theory  as  does  the  adjustment  for  normally 
distributed  area  targets.  Probabilities  of  damage  calculated 


using  these  adjustments  will  usually  be  within  2 or  3 percentage 
points  of  the  rigorously  calculated  results. 

These  CEP  adjustments  can  be  easily  incorporated  into 
the  computer  program  of  Appendix  E for  rapid  calculations  of 
probabilities  of  damage. 

C.  Single  Shot  Probability 

Single  shot  probability,  SSP,  measures  the  relative 
frequency  with  which  a weapon  can  be  placed  within  a designated 
area  assuming  that  with  a large  number  of  independent  trials 
there  would  be  a circular  normal  distribution  of  impact  points . 
For  a circular  normal  distribution,  the  probability  that  an 


impact  point  falls  in  an  infinitesimal  area  dA  a distance 


from  the  DGZ  is 


I 


f 


dA.  As  discussed  in  section  V, 


2rr<rl  ® 

tne  standard  deviation  of  this  distribution,  C T,  is  equal  to 
.84932  CEP.  The  SSP  is  found  by  integrating  this  density 

c 

function  over  the  designated  target  area,  A. 


If  the  target  area  is  circular,  it  is  easiest  to 
perform  this  integration  using  circular  cylindrical  coordi- 
nates (r,9)  centered  on  the  target.  The  general  case  where 
the  aim  point  is  offset  a distance  x from  the  center  of  a 
target  of  radius  TR  is  shown  in  figure  6. 
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the  relationship  between  j?  and  r is 

f 


z = r*  + xz-  Zra  c<*  © 


The  SSP  is  then 


SSP  = 


2TT<r 


Zir  tr  2.*"*  cos  © 

Z<r* 

a.  \ l G 


ZTI  T 

n 


rdrd© 


= e 


T /jno  I ,'Z*  ^Ccse 
where  _L0V.<r"-;  — ZtFJ  ® d0  is  the  zeroth  order  modified 

Bessel  function  of  the  first  kind  discussed  earlier  in  this 

section.  -• 

If  the  aim  point  is  located  at  the  center  of  the 


circular  target  (i.e.,  x=o)  then  the  above  equation  reduces  to 

SS  P = ;p.  f\  erV<lr 

o 

a 


- I - e-s(TR/r) 
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In  section  V the  CEP  is  defined  such  that 
This  definition  can  be  used  to  eliminate 
expression  for  SSP  to  obtain 


_L  = I - e 
z 1 . 

CT  from  the 


above 


5S  P = I-.5 


(t  r/cep) 
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D.  Equivalent  Target  Area 

For  certain  special  classes  of  targets  such  as  bridges, 
dams,  locks,  runways,  etc.,  a specified  degree  of  damage  to  some 
part  of  the  target  satisfies  the  damage  objective.  For  example, 
for  a bridge  the  collapse  of  one  span  is  usually  the  damage 
objective.  The  exact  determination  of  the  probability  of 

c c , 

damaging  such  targets  is  possible  but  quite  laborious.7 
Approximate  answers  can  be  obtained  with  much  less  effort  by 
employment  of  the  equivalent  target  area  (ETA)  approximation. 

The  ETA  is  defined  as  an  area  such  that  the  probability  of 
placing  the  GZ  in  the  area  is  equal  to  the  probability  of  doing 
the  desired  level  of  damage  to  the  target. 

For  a rectangular  target  the  ETA  is  approximated  by 
adding  marginal  strips  around  each  edge  of  the  target  of  width 
equal  to  the  weapon  radius  for  that  aspect  of  the  target.  Then 
tne  probability  of  damage  to  a rectangle  having  a length  Z and 
width  w is  approximately  the  probability  of  placing  a weapon  in 


-7"Target  Analysis  for  Atomic  Weapons",  p.  5-1.  * ■ 
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a rectangular  area  of  length  (£.+2WR^)  and  width  of  (w+2WRw) 
where  WR^  is  the  weapon  radius  associated  with  the  length  VN 
and  WRw  is  the  weapon  radius  associated  with  the  width  VN. 

This  approximation  is  improved  if  the  actual  CEP  of 
weapon  delivery  is  replaced  by  an  adjusted  CEP  for  each  aspect 
of  the  target  whi*ch  compensates  for  the  variation  of  damage 
probability  with  distance.  The  width  adjusted  CEP,  CEPW , is 


given  by  CEPW  * 1.1774  aw  = /CEP2  + (1.1774  ad)2  WR^2 , and  the 
length  adjusted  CEP,  CEP^  , is  given  by  CEP^  = 1.1774  a ^ = 

^CEP2  + (1.1774  ad)2  WR2  where  cd  is  the  damage  sigma  having 
a value  between  0.10, and  0.50,  and  ow  and  a ^ are  the  adjusted 
standard  deviations  of  the  adjusted  delivery  function. 

If  the  DGZ  is  located  at  the  origin  of  the  x,  y 
coordinate  system  of  figure  7,  the  probability  that  a weapon 
falls  in  the  semi-infinite  plane  from  x = a to  x = oo  is 
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Figure  7 
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t = Holding  "x"  constant  and  integrating  with  respect  "y”  yields 


Let  t - x/0|,  dx  = o^dt.  Therefore, 


where  |a| is  the  absolute  value  of  a,  i.e.  |-5|  = 5. 

The  probability  that  the  weapon  falls  in  the  semi- 
infinite  strip  of  width  (b-a)  is  (Pa-Pb)  , where  b-a  = l + 2WR£ 
(see  figure  7).  The  probability  that  the  weapon  falls  in  the 
semi-infinite  strip  of  width  (d-c)  is  (Pc~Pd) , where  d-c  = 
w + 2WRw.  Therefore,  the  probability  that  the  weapon  falls  in 
the  rectangle  defined  by  x = a to  b , y = c to  d is  P , where 
p “ <pa--V  <pc-pd>- 


I 

£ i 


This  result  is  easily  programmed  using  the  series  expansion 
for  the  error  function  from  Section  IV.  A. 

The  ETA  for  a circular  target  is  constructed  by  adding 
a ring  of  width  equal  to  WR  around  the  target.  The  probability 
of  damage  is  then  calculated  using  the  method  of  part  C,  Single 
Shot  Probability , where  TR  is  taken  to  be  the  target  radius, 
and  the  CEP  is  replaced  by  an  adjusted  CEP  (CEPa)  where 
CEPa  = /cep2  + (1.1774  ad)2  WR2  . 


VII.  VULNERABILITY  NUMBER  (VN)  CODING  SYSTEM 

A.  Vulnerability  Numbers 

In  the  VN  system,  a target's  susceptibility  to  blast 
damage  is  indicated  by  a combination  of  numbers  and  letters. 

The  vulnerability  number  (VN)  consists  of  a two-digit  number 
reflecting  the  target  hardness  relative  to  a specified  damage 
level,  a letter  indicating  predominant  sensitivity  to  over- 
pressure (P)  or  dynamic  pressure  (Q) , and  a K factor.  The  two 
digit  numerical  value  scale  of  the  VN  is  an  arbitrary  classi- 
fication describing  a target's  hardness.  It  is  a linear  function 
of  the  logarithm  of  the  peak  pressure  from  a 20  KT  weapon  that 
would  have  a 50%  probability  of  damaging  a randomly  oriented 
target  to  the  desired  level.  The  base  yield  was  chosen  to  be 
20  KT  instead  of  the  more  convenient  1 KT  because  the  original 
system  was  developed  from  the  Hiroshima-Nagasaki  data  assuming 
that  tne  yields  of  the  Hiroshima  and  Nagasaki  weapons  were  20 
KT.  The  appropriate  damage  sigma  for  P targets  unless  other- 
wise specified  is  (Jj  = .20.  The  appropriate  damage  sigma  for 
Q targets  unless  otherwise  specified  is  0^  = .30.  The  K factor 
allows  for  hardness  adjustments  to  be  made  to  account  for  the 

C ; 

effects  of  variations  in  blast  wave  duration  due  to  different 
weapon  yields.  Each  VN  must  also  have  a specified  damage- 
level  criterion,  such  as  "collapse, " "24-hour  recovery  time," 
"severe  damage  to  contents,"  "moderate  structural  damage, 


II 


etc. 


The  completely  arbitrary  coding  relationship  for  p 
type  targets  is 


PVN 

I.I2IC*  Cl.  2) 

where  pj#  is  the  peak  overpressure  from  a 20  KT  weapon  that  has 
a 50%  probability  of  achieving  the  desired  level  of  damage  to 
the  target.  This  relationship  was  established  so  that  a p 
of  10  psi  corresponds  to  a PVN  of  12.  With  each  integer  in- 
crease in  the  PVN  the  associated  p50  increases  20%. 

This  ceding  relationship  may  be  inverted  to  obtain 


or 

PVN  - I2..L3  — .H 

Since  the  peak  overpressure  at  a given  range  is 
uncertain  to  roughly  * 20%,  this  coding  relationship  insures 
that  P type  target  hardnesses  are  not  specified  more  precisely 
than  justified  by  the  pressure-range  data.  This  scale  con- 
veniently allows  for  the  complete  pressure  range  of  interest 
to  be  coded  by  a two-digit  number. 

The  dynamic  pressure  coding  scale  was  chosen  using 
the  approximate  form  of  the  Rankine-Hugoniot  equation, 
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q = .023  p2.  The  scale  was  defined  So  that  the  dynamic  pressure 
required  for  a 50%  probability  of  damage,  q50,  for  the  VN  of 
interest  is  equal  to  .023  P502  where  ps0  is  from  the  numerically 
equal  P VN.  Therefore,  the  VN's  are  "tied”  at  the  50%  proba- 
bility but  not  at  other  probabilities.  Thus, 


qso  = .023  (pso2) 

= .023  (1.1216*(1.2)QVN) 2 
q 5 0 = .02893*(1.44)QVN 

The  above  relationship  establishes  a 44%  difference  :n  peak 
dynamic  pressures  between  adjacent  VN's.  This  can  be  inverted 
to  give 

Q VfJ  h°  ~ 

; ‘V  ■ ■ ■;  , ■ : ?*v  ,‘c.  ■ ;v',  / 

or 

QVM  = £.31  Ay.  Z 

This  relationship  insures  that  the  complete  dynamic 
pressure  range  of  interest  is  covered.  Since  the  dynamic 
pressure  at  a given  range  is  only  known  to  within  about  + 40%, 

Q type  target  hardnesses  are  also  not  specified  more  accurately 
than  justified  by  the  pressure  range  data.  • ” 1 : 
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Care  should  be  taken  not  to  use  the  equation  q = 

.023  p2  to  calculate  the  peak  dynamic  pressure  corresponding 

to  a given  peak  overpressure.  For  overpressures  below  about 

10  psi  the  equation  q = .023  p*  is  a good  approximation  to  the 

5 -fl? 

Rankme-Hugomot  relation  c^=  — where  p0  is  the  ambient 

atmospheric  pressure.  This  Rankine-Hugoniot  relation  was 
derived8  assuming  an  ideal  shock  front.  It  only  fits  the 
available  experiment  data  fairly  well  for  zero  heights-of-burst 
(HOB)  . The  correct  determination  of  the  q given  p or  vice 
versa  for  a given  HOB  must  be  through  the  horizontal  ground 
range  using  pressure-range-HOB  curves  such  as  figures  1-4  and 
1-9  of  AP  5 50-1-2-69 -INT . 

B.  K Factor 

As  previously  mentioned,  the  blast  wave  duration 
varies  with  weapon  yield.  The  increased  blast  duration  associ- 
ated with  larger  yields  may  cause  targets  to  fail  at  lower 
pressure  levels,  while  at  small  yields  the  reduced  blast  dura- 
tion may  necessitate  higher  pressures  for  target  failure.  To 
account  for  this  yield  dependence/  the  PV  system  uses  K-f actors 
for  both  P and  Q targets . The  K-f  actor  is  an  integer  from  0 to 
9 which  adjusts  the  base  VN  to  reflect  the  sensitivity  of  the 


8Samuel  Glasstone,  "The  Effects  of  Nuclear  Weapons",  Air  Force 
Pamphlet  No.  136-1-3,  Department  of  the  Air  Force,  April  1962. 
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target  to  the  different  pressure-time  pulse  shapes  for  yields 
other  than  20  KT.  A K factor  of  0 indicates  a target  that  is 
not  sensitive  to  blast  wave  duration  and  can  be  expected  to 
fail  at  the  same  pressure  regardless  of  weapon  yield.  A K 
factor  of  9 indicates  a target  that  is  very  sensitive  to  blast 
wave  duration  and  can  be  expected  to  fail  at  quite  different 
pressures  at  various  yields . 

The  adjustment  factor  R is  the  ratio  of  the  pressure 
(either  overpressure,  p(Y),  or  dynamic  pressure,  q(Y),  required 
for  a 50%  probability  of  damage  at  yield  Y to  the  pressure 
required  at  20  KT  (p(20)  or  q(20)).  The  K factor  is  related  to 
the  adjustment  factor,  R,  in  the  following  manner^ 


o _ I K ( . \ _ ffilL 

R - I ~ io  td  j ~ 


or 


±iXL 


where 


tdo  = positive  phase  blast  wave  duration  for  20  KT 
td  = positive  phase  blast  wave  duration  for  yield  Y 

* Y 

t,  s .45  * for  overpressure 

. . . « ' 

* — for  dynamic  pressure  . 


td  s .105 


f 


"Development  of  the  K-f actor  in  the  VN  System,"  PV-105-61, 
Air  Force  Intelligence  Center,  Assistant  Chief  of  Staff, 
Intelligence,  25  January  1961. 


Therefore,  for  P type  targets. 


K- '-***(-*& 


K o'/*-  ( 20^* 


For  Q type  targets. 


R=I-^+-^(-|S)/3(t)3=  '°+  • 


These  equations  are  solved  for  R for  a given  K factor  and  yield 
by  iteration  to  the  desired  accuracy  using  a first  guess  of  2 
for  PVN's  or  a first  guess  of  3 for  QVN's. 

The  adjustment  factor  R is  used  to  determine  the 
adjusted  VN  (VNa)  using'  the  PVN  and  QVN  coding  relationship  as 
shown  below.  por  P Vtl’s,  VNa  = VN+A,  where 


A = j-?,g.w 

log  (1.2) 

« 12.63  logl0  R 
* 5.485  loge  R 


For  Q VN's, 


log (R) 
log (1.44) 


* 6.315  logi#  R 
« 2.742  loge  R 


C.  VN-Probability  Relations 

It  is  often  necessary  to  know  not  only  pso  for  a 
given  VN  and  yield,  but  also’ the  pressures  for  other  proba- 
bilities of  damage.  The  method  used  to  obtain  the  pressure 
p,  for  a probability  of  damage  a%  for.  the  adjusted  VN  of  v2  is 

a "" 

discussed  below. 

For  P type  targets , the  VN  coding  relationship  gives 

VZ 


' P50=  1-121^  0-2) 


The  analysis  of  the  Hiroshima-Nagasaki  and  Nevada 
test  data1®  resulted  in  the  adoption  of  the  following  relation- 
ship between  the  overpressure  required  to  damage  a structure 
to  a given  level  and  the  pressure  which  gives  a 50%  probability 
of  damage,  p5Q. 


.217  b 


-p.  = 'Pso  e 

' | 'C  ' t ■ -='\r  _ ; ' ' . c , 

The  probability  a is  given  by 

,b 


a - 


1 

-fzrx 


{ 


<4x 


where  b is  the  probability  a expressed  as  probits-5.  (A  probit 
has  the  magnitude  of  the  standard  deviation.  Minus  5 probits 
is  defined  as  a = 0%,  5 probits  is  50%,  and  10  probits 


*■» 


1®"A  Classification  of  Structures  Based  on  Vulnerability  to 
Blast  from  Atomic  Bombs",'  - l . . . ' ■ <• 

''  V ^ ‘ , . ' "c  40  . \ , c --  o ° .C 


I 


. corresponds  to  a = 100%.) 

» - - Similarily,  for  dynamic  pressure, 

1 vi 

=.02813  ( 1.44) 

I.04Z  b 

%fl.  ~ so  e 

;■  ;•  • 

In  addition  to  coding  P50's  and  q^'s,  adjusted  VN's 
, code  the  scaled  weapon  radii,  SWR11,  to  be  used  in  calculating 

Pd(r),  the  distance  damage  function.  SWR  for  each  VN  are 
£ calculated  using  the  pressure  - probability  functions  given 

\ above  by  combining  them  with  pressure-range-HOB  data.  SWR's 

calculated  in  this  manner  are  listed  in  AP-550-1-2-69-INT  for 

\ \ ■ 

' PVN's  and  QVN's  at  various  scaled  heights  of  burst,  SHOB.  To 

, facilitate  computer  applications,  these  calculated  SWR's  have' 

i been  approximated  with  polynominals  of  the  form  SWR  = f(v2) 

I for  various  scaled  heights-of-burst . The  polynomials  have  the 

( - 
I 

! form 

(n  • 

E k:  ( vzy 

y°  i 

! where  n varies  from  2 to  7 depending  on  the  SHOB.  The  con- 

, * stants  and  a further  description  of  the  curve  fits  are 

1 3 

i -r-r 

< < ■‘■■‘■Scaled  weapon  radii  and  scaled  height-of-burst  have  dimensions 

f t/KT^i * 3 . 
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Numerical  Calculation  of  Momenta 


The  moments  presented  in  Section  III. A.  can  be  calculated 
numerically  if  the  damage  function  Pjj(r)  is  known  at  the  ranges 
r^  either  from  experimental  data  or  from  the  analytic  approxi- 
mation to  P^fr)  from  Section  II. B. 

Consider  the  contribution  to  the  first  moment  from  the 
small  section  of  P(j(r)  from  r^_^  to  r^  shown  in  figure  A.l. 


r 


Figure  A.l. 

Over  the  small  range  from  r^_^  to  r^ , accurate  estimates  for 
P<j(r)  can  be  obtained  by  linear  interpolation  between  the 
values  of  P^  at  r^_i  and  at  r^.  Therefore, 

Pd(r)  » P,(n)-(Pjr,)- 

The  contribution  to  from  this  section  of  the  damage 
function  is  p 

<r>;  = { -(PM - pd fn-OXo^r)]  dr 

ri-l 

‘ 1 • . ’ 43  Appendix  A 


<rV*T(e<ra+fc<rJXn-  *i-.)  ■ ' 

Since  the  first  moment  is  the  summation  of  all  possible  terms 
like  that  above, 


where  rQ  = 0. 

Similarily,  the  second  moment  (r*}  is  approximated  by 


<r*>=!  i[R,('-.Xn-n-,)(ri.,  + f (r;-r;.,))  + 


(P/O-  P,CO)(!1ykl-)(ri.,+  y(r1- 


This  form  may  be  rewritten  as 


where  again  rQ  = o.  ‘ 

As  explained  in  Section  II. A.,  the  weapon  radii  (WRa^r *•>*,) 
presented  in  Part  I of  AP-550-1-2-69-INT  are  calculated  numeri- 

c •'  c , - : , . «,  ■ , 

cally  using  this  formulation.  The  damage  sigmas  ( (%  = y | - ) 

calculated  numerically  for  P type  targets  range  from  .10  to  .20 


■«vr»  *-  -ttHtKfv’nt.,- ~ 


while  for  Q type  taryets  the  damage  sigmas  range  from  .20  to 

I .30.  For  probability  calculations,  the  distance  damage  function 

is  approximated  bv  the  analytic  function  of  Section  XI. B.  The 
parameter  WR  of  this  function  is  taken  to  be  the  numerically 
calculated  WR,  while  the  parameter  is  assumed  to  be  .20  for 
P targets  or  .30  for  Q targets.  This  analytic  approximation 
to  Pjj(r)  will  always  have  P<j(o)  = 1.  However,  experimental 
data  may  not  have  Pj(o)  = 1 due  to  the  hardness  of  the  target 

I;  and  the  particular  yield  and  HOB  chosen.  In  such  cases, 

j 

t accurate  probabilities  of  damage  can  not  be  obtained  using 

i 

j the  analytic  approximation  of  P^fr). 
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Derivation  of  Combined  Effects  Formulas;  Section  IV 


The  derivation  of  Equation  1 from  Section  IV  is  as  follows: 
From  Section  IV, 

WR*  = \A/R*  t WR*  - 2 fr  Pd( r)  Pd  (r)  dr 

■ii 

•• 

Let  I = Z (r  Pd  (r)  p.  (r)  dr 
Using  integration  by  parts,  let 

u “ 6/0  Rix(r)  and  dv  = 2rdr 

then  du  = Pd  (r)  d Pj  (r)  + p,  (r)  d p,  <» 

» 31  X i 

and  v=r2  . ‘ 

I then  becomes 

#9  CO  0Q 

I = -{r’PjMdPjW-Jr^PjWdPj.Cr) 


. = - J r2(J  (r)d  P,  (r)  - /V*  Pd  (rl  d PjM  (EQU.  B-l) 

° 1 ' * • ' * 

Since  the  terms  are  symmetrical  in  Equation  B-l,  it  is  necessary 
only  to  deal  with  one  and  apply  the  findings  to  both  integrals. 


Let  3L=  - | d 

o' 

From  Section  II. B.  we  know  that  fj(r)=  j'  (p(lj)^lj 


2,(r) 


where  (p(y)  - G , 
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Since  d(isn(r))  - and  r*  = Q 2£*(r'  , the  integration  of 

Equation  B-2  may  be  made  with  respect  to  ln(r)  rather  than  r. 

lim  ln(r)  = oo  lim  ln(r)  = -oo 

r-*oo  r-*o 


In  (r ) 


Combine  the  exponential  terms  and  complete  the  square  to  obtain: 

r / \ / , „ \i 


ii  = e 


/ ( £ 4>(y)  ex  P f ^ r " ^ (Wl *z )}) 


d 

-iTrf  r- 


If  we  let  y 


X — ty\  r 

f9. 


and  dy  = 


, then 


47 


ii  = e 


ln(wR 


tn(\ug,£?‘ ) 


H 

S’ 


dfcir 

j3a 


Interchange  the  integral  siyr.s , combine  exponential  terms  and 

complete  the  square  to  obtain 

4*(tuR,<f  ) 


* 

f where  B = 


C = 


x&  + im  (vjRze^x ) 

'f  Z5 


and 


P m P**f 3 


x 

X 


Holding  the  variable  "x"  constant  and  integrating  with  respect 
to  ln(r)  yields  a value  of  -jg-  for  the  inner  integral . There- 
fore , 

• „ ’ -V  ••  ‘ 


J 'YarF  exp["i'/3i^(c"SZ)] 


Simplifying  the  exponential  under  the  integral  sign  this  becomes 
• V . U(urI)  ' z 

~~P  : J ' YXF  eX?["i^(x"^(WR*e^)  Jdx 


II  = 
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a* 

. , _ X-  ^y\(wRxerx) 

Let  j?  = , then 


ii  = e 


4n  (WR*) 


( 4(f)  df 


where  cL,=  C ^ ) = 2^  ^r- WRie^x^  • We  then  obtain 

II  = WR  \ 

. oo 

Employing  the  same  procedure,  the  other  integral  in  Equation  B-l 
is : 


WRi  ( 4>($)  d§ 


where  ^ = j jU  ^ ^ ) s Z1  ( r ’ WR  e )* 

Substituting  into  the  equation  for  WR?  2 from  Section  IV, 

Otj,  *■, 

wr^2  = wrT+w rI  - WR,1  C '</>(§)<!$  - WRZ  ( 

WE-’  - WR,Z(|-  + WRZ  ( | - 


where 


and 


OU 

j = B(z'(  WRzePl)) 

\<p(e,)di,  =B( z;(WR,e^))  . 
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This  result  is  equation  1 of  Section  IV. 

The  mean  for  combined  weapon  effects  as  presented  in 
Section  IV  is 

OO 

< r'*>  = f ( + pJiCr)  - P4(. r)  Pj/O)  dr 

OO 

-<>-,>+<  f Pj/O  P^Cr)  dr 

OO 

Let  1 = / P,'(r)  PjJr)  dr. 

Using  integration  by  parts,  let 

u = fj'(r)  Pj%(r)  and  dv-dr. 

Tnen  du  ■=  Pj/r)  d + Pj/r)  d Pd/r> 

and  v=r. 

Our  integral  is  then  , - ..  ; 

«o  pp 

i = - if  r f^COdPj^r)  - £r  pdjr)  d . 

Since  the  integrals  are'  symmetric , it  is  sufficient  to 
examine  only  one  integral  and  apply  the  results  to  both. 

•90  c 

Let  II  = - f r PitCr)'AP^ Cr) 


where 


50 


We  also  know  that 


dBuCO  dr  _ ~eXPL-ip-(/n(wi 

dr  ~ 'iti.Tf 

Therefore,  II  becomes 

.A/  w r i > 


r-jir(^CwR2e^)-A  r)A] 


ii  -J z)-/n  r)  J jd[ 

/*i  r 

Changing  the  variable  "r"  to  ln(r)  such  that  r=  e an<3 

d £nr  = -^r-  , this  becomes 


» -P((i  zn.y*. 


Combining  exponential  terms  and  completing  the  square, 


rexpL“zf 


r — An  W(?z 


d ioa 


X — An  r Jy 

Letting  y = then  dv  =—3 — and  when  y=z  (r)  then 

j/3'  r»  » 

x = In  (WR(  e^'  ) . Therefore, 


t(wR,e  ) 


II  = <(V. 


Combining  exponential  terms  and  completing  the  square 

t t _ / n v f f(  ( ex  p C-  -irCfe-Yi (c-  BZ)J  , Y ex  p[“  (fcir-B)_. 


=<r# 


ex  P C-  &)±  jY  ex gkjAgk 

VZJr  /s,  71  nrarr  jS* 


where 


B 


Xftx  + In  Wl?2 

P - 

- + Is, 2(  /*  wgaf 

r 


- /» i 4 ^ 


and  /S-=|8-+/3- 

Holding  "x"  constant  and  integrating  with  respect  to  ln(r) 

Ui(wa,e^  ) 

ii  - <ra>  ( S'tl-frl*-**'”**)*}  dx 

J V^T 


Letting 
II  = 


- if -*$- 


then 


<ra>{-vfeef/^df  where  *•  = yM*twtr ) 


Therefore,  the  first  moment  for  combined  effects  is 


*•  /* . 

<r-*>  “<r«>J  -yfr  e 


-S% 


= <r.>l 

oc 

^eS/2d5  + <ra) 

t 

*» 

H 

II 

/\ 

b(*;(wr4))j 

where 

B ( 2 - (rj)  =£^6  dt 

This  result  is  equation  2 cf  Section  IV. 


Calculation  of  Weapon  Radii  for  Combined  Effects 


As  explained  in  Section  IV  when  two  or  more  weapon  effects 
significantly  contribute  toward  damaging  a target,  a combined 
WR  and  should  be  calculated  and  the  P.D.  'computation  then 
based  on  this  combined  WR  and  a^.  This  procedure  is  most  often 
used  in  estimating  personnel  casualties/fatalities  when  radi- 
ation and  blast  have  comparable  ranges  of  effect. 

The  example  program  takes  two  WR's  and  their  associated 
damage  sigmas  and  calculates  the  WR  and  damage  sigma  of  the 
combined  effects. 

WRi  2 , the  combined  WR,  <ria>,  the  combined  mean  radius, 
and  od  j , the  resultant  damage  sigma  of  two  damaging  effects, 
may  be  calculated  using  the  equations  of  Section  IV. 

wr5*  = WR*^  e3^  + WR*  fy==  e 

22(wR,ePl  ) *,(wR3tePl) 

<ra}=  <r,>  ( 7=^  e +-<r*>  C^=.  e ^ du 


adi2  ~ I — 


<r,*y 

VJR~ 


I z 


Pm 


wne 


I t / MiA.  \ 
re  Z;(d)=  Yt  ( d ) 
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In  the  example  program: 

f -s*. 

f i = ] e 

f2  = < e~^2dy 

*2. 

and  fna(v)  = erf  (v) 

t 

t is  an  indexing  variable,  if  t=0  then 
zl  = 2X (w2) 
z2  = z2 (wl) 
fl  = fl(z, (w2) ) 
f 2 = f2 (z2 (wl) ) 

and  r = <r  , 2>  = vl*exp(-^2/2)  *f 2+w2*exp (-3z/2)  *fl 
since  r^  = WR^*exp  (-0  z/2) 

After  r is  calculated  then: 
zl  ■ zt (w2*exp (0Z) ) 
z2  = z 2 (wl*exp  (f3z) ) 

t is  then  set  to  1 and:  ~ fA.  <■/  i 

fl  = f 1 (z  ^ (w2exp ( 8 2 ) ) ) 
f 2 = f 2 (z  2 (wlexp ( B 2 ) ) ) 

and  w = WRj  2 = ywlz*f2  + w2z*fl 

After  r and  w have  been  calculated,  tHe  combined  damage 
sigma  may  be  derived:  r ’ 


05  Tt.m  COMPUTES  THE  COMBINED  VP  AKr  SIC-MA  : Given  2 wr  and  s 
10  data  .07052  50784,  .04228  20125,  .00927  05272 
20  data  .00015  20143,  .00027  65f??2,  .00004  .'SO 638 

30  read  a,b,c,d,e,f 

40  def  fna(x)  = 1 - l/(l+a*x+b»x ?.+c*x'Z+&*x  '4+e*x  5+f*x  6)'16 
50  print’input  vr(l),s(l),wr(2),s(2)“, 

60  input  wl,sl,w2,s2 
65  let  t = 0 

70  let  bl  = sqr{-  log(l-sl^2}) 

80  let  1)2  = sqr(-  1 04? ( 1 — s 2^2 ) ) 

90  let  1x3  = sqr(bl  2+b2 *2) 

100  let  z 2 * (l/b3)*lo?(w2*exp(-b2  2)/wl) 

110  let  zl  = (l/b3)*lo^(wl»exp(-'bl'2)/v2) 

120  let  fl  = ,5-.5#zl/abs(zl)*fna(abs(zl)/1.41421) 

130  let  f2  = .5  - ,5*z2/abs(z2)*fna(abs(z2)/1.41421) 

135  if  t > 0 then  180 

140  let  r * vl#exp(-bl~2/2)*f2  + w2»exp(-h2"2/2)*fl 
150  print  ”<r>  = ’’  int(r+.5) 

160  let  z2  = ( l/b3)#l04( v2*exp(-b2  2)/wl/exp(bl  2) ) 

170  let  zl  = (l/b3)*l04(wl*exp(-bl  2)/w2/exp(b2  2)) 

173  let  t = t + 1 
175  r 0 to  120 

180  let  v = vl"2*f2  + w2*’2*fl 
190  let  v = sqr(w) 

200  print  "VR  = "int( w+.5) 

210  let  s = sqrtl-r  2/v*'2) 

220  print  "si^ma  *5  ”int(l00*s+.5)/l00 
230  print 
240  print 
250  ro  to  50 


N K \ 

\ f ! 

\ . . REAIJY 

.■  \ * 

•run 


input  wr(l),s(l),wr(2),s(2)  ? 

3500,. 3,3000,. 4 

<r>  = 3709 

TO  = 3P73 

sif'fita  = .29 


Input  vr(l)  ,s(l)  ,vr(2),s(2)  ? 

10500,. 4, 8000,. a 
<r>  « 10644 

TO  = 1137ft 

sl4»a  = .35 


input  wr(l), fi(l), wr{2), s(2)  ? 

5000,. 4, 3500,.  *■ 

<r>  * 4953 
TO  *=  5306 
slfwa  - .36 


Input  vr(l),s(l),*r(?.),fi(2)  ? 

j 1200,. 3,1000,. 3 

<r>  « 1245 
f ■'  TO  - 1289 

5 sl^ma  - .26 


t 


74/07/02  tibh  1 ->• 
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i. 

I \ 

Weapon  Radii  Determination  for  P and  0 
Targets  Including  K-Factor  Adjustments 

Scaled  weapon  radii  (SWR)  for  PVN's  and  QVN's  numerically 
calculated  at  various  heights  of  burst  as  explained  in  Appendix 
A may  be  approximated  by  polynominal  curve  fits  of  the  form 

log  (SWR)  = J L (V£)i 

.j*o 

where  v2  = adjusted  VN 

and  n varies  from  2 to  7 depending  on  the 
scaled  height  of  burst 

The  constants  kj  are  calculated  for  SWR's  associated  with 
PVN's  and  QVN's  at  each  100  f t/kT  ^3  increment  in  SHOB  from 
, surface  to  900  ft/kT^3.  Fits  at  other  HOB's  are  not  provided 

because  SWR's  were  calculated  only  at  100  ft/kT^3  increments 
in  scaled  HOB.  This  increment  supplies  sufficient  accuracy 
when  interpolating  for  a non-tabulated  HOB.  Fits  are  not 
provided  for  higher  HOB's  because  higher  HOB's  would  be  above 
tne  "knee"  in  the  HOB  curves  for  pressures  of  interest. 

For  PVN's  at  scaled  heights  of  zero  ft.,  100  ft,  800  ft. 
and  900  ft.,  two  equations  per  SHOB  are  required.  These  are 
tnird  and  fourth  degree  polynominals  with  one  fitting  VN ' s 
<7.5  and  the  other  fitting  VN's  >7.5.  However,  at  all  other 
overpressure  SHOB ' s and  at  each  dynamic  pressure  SHOB,  one 
equation  gives  an  accurate  fit.  These  polynominals  are  of  the 
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sixth  or  seventh  degree. 

One  advantage  to  this  polynominal  fit  is  that  it  gives 
a better  interpolation  for  fractional  adjusted  VN's  than  can 
be  achieved  by  interpolating  between  SWR's  as  given  in  a list . 

A disadvantage  is  that  optimum  SHOB’s  must  be  chosen 
even  to  100  ft.  If  the  true  optimum  SHOB  were  150  ft.,  the 
program  would  provide  for  calculating  the  WR  at  100  ft.  and 
at  200  ft.  and  then  interpolate  between  these  two  to  arrive 
at  a WR  for  150  ft.  Thus  if  100  ft.  and  200  ft.  were  less 
than  optimum,  the  WR  computed  at  150  ft.  will  also  be  less 
than  optimum.  Therefore,  to  obtain  the  largest  WR  possible 
using  the  given  curve  fits,  select  the  SHOB  even  to  100  ft. 
wiiich  gives  the  largest  WR.  This  calculated  WR  will  generally 
not  differ  from  the  actual  optimum  WR  by  more  than  5 ft/kT1^. 

The  following  program  is  designed  so  that  VN,  type  target, 
K factor,  appropriate  damage  sigma,  yield,  and  height  of  burst 
are  input  for  each  case.  The  appropriate  damage  sigma,  s,  is 
input  to  enable  WR’s  to  be  computed  for  targets  which  do  not 
have  SIGMA-20  associated  with  their  PVN  or  SIGMA-30  associated 
with  their  QVN. 

, The  VN  adjustment  is  performed  using  the  following 
variable  definitions. 

vl  = base  VN 
k = k factor 


Y = yield,  in  kilotons 

v = adjustment  to  the  base  VN  = A*log(R) 
v2  = vl+v  = the  adjusted  VN 
The  adjustment  to  the  VN  is  derived  by  solving  the 
equation 


R*  I - (is)  + 3 


(EQU.  D-l ) 


where  e = 1/2  for  PVN 
e = 1/3  for  QVN 

The  solution  for  R may  be  achieved  by  an  iterative  process. 
Modify  Equation  B-l  to  read 


n = I 


(EQU.  B-2 ) 


Initially  set  PI  = 2 if  PVN 
n = 3 if  QVN 

Solve  for  r2.  If  |r2-rl|  > .001  then  let  rl  = r2  and  again 
solve  Equation  B-2  for  r2 . Continue  this  process  until 
)r2-rl|  < .001 

The  adjustment,  v is  equal  to  A*loge(r2) 
where  A = 5.485  for  PVN 
A = 2.742  for  QVN 
The  adjusted  VN  is  then 
v2  = vl+v 
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The  WR  is  calculated  using  the  following  variable 
definitions. 

h = actual  height  of  burst 

lA 

hi  = scaled  height  of  burst  = h,/Yr 

h(l)  = int (hl/100) *100 , i.e.  the  nearest  even 
100  less  than  hi 

h(2)  = h(l)+100,  i.e.  the  nearest  even  100  greater 
than  hi 

w(i)  = scaled  weapon  radius  at  h(i) 
w = SWR  at  hi 


wl  = WR  scaled  to  correct  yield  and  damage  sigma 

s = appropriate  damage  sigma,  .2,  .3,  .4,  .5 

k0,...k7  = appropriate  constants  for  the  polynominal  curve 
fit;  stored  for  each  100  SHOB  up  to  900  ft. 

A scaled  weapon  radius  is  computed  at  h(l)  and  at  h(2). 

If  hi  = h(l)  then  a computation  is  made  only  at  h(l) . 

m • 

g(i>  - ko  +■  IT  k-  (vzy 

j"  * 

w (i)  = ex  p (g(i)) 


The  SWR  for  hi  is  derived  by  linearly  interpolating  between 
w(l)  and  w(2)  . = 

, w=w(l)+e*(w(2)-w(l)) 

where  e = (hi  - h(l))/100 

W is  then  scaled  to  the  correct  yield  and  damage  sigma 
to  give  the  actual  WR.  SWR's  for  PVN's  are  based  on  a damage 
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\ ‘ • 

| sigma  of  .20  and  SWR’s  for  QVN's  are  based  on  a damage  sigma 

X 

of  .30.  If  the  appropriate  damage,  s,  differs  from  this 
generality,  the  suitable  adjustment  is  based  on 
j WR  = r 50/ (1-s  2 ) 

si  - 1.04  = 1/(1-. 22)  for  PVN's 
■ . si  = 1.10  = 1/(1-. 32)  for  QVN's 

Thus,  w * Y^/sl  = r50  = range  such  that  Pd(rso)  = .5 
al  = l/(l-s2) 

wl  - r50  * al  = w * yVz  * al/sl 


♦list  • 

OF  rem  KVRCAL- Calculates  Weapon  Fadii 

10  nrinfllOTT  VK,T,K,?,yieU,Kob" 

20  input  vl,t?,k,$,y,h 
on  if  t?  = '>'■  then  PO 

40  let  e = (i/0>) 

50  let  si.  = l/(l-.3*2) 

00  let  a = 2.742 

ft  le  t n - 

PP  ro  to  ISO 

PO  let  e = .5 

1CP  let  si  * 1/(1-./  0) 

110  let  a = 5.405 
1?0  let  rl  * ? 

IOC  let  r?  = l-(k/in)+(k/10)*(?r./yV'(l/?5)*(rl)  f 

V-f  if  abf  (rt-rl)<  .CfJ  then  l^n 

IfC  let  rl  = t2 

150  ro  to  100 

170  let  v = a*lor(r2) 

1°0  let  v?  s vi+v 
IPO  let  hi  = h/y"(l/3) 

200  le  t 1.(1)  «int(hl/icr)*100‘ 

210  for  i = 1 to  2 
220  if  tf  - * p"  tl  en  ?PC 

77 1 ro  to  1*30 

272  let  r(i)  = k.O  +kl*v2+k2*v2'2  +k?*v2  7 ♦ K4*v2  4 
274  ro  to  250 

c?40  let  r(i)  * k 0+kI*v?.  +k?*v?"  2+k.3*v2  3+k4*v2  4+k5*v2  5+k.e*v2  H+k7*t2  V 
7*0  lr  t v(l)  = e*p(r(i)) 

2fO  if  atsfhl-h(l)  ;<  .01  then  2R0 
< 270  ro  to  7% 

290  let  i =2  ' - 

2P0  let  e = 0 

30C  next  i 

710  let  al  = 1/U-s  2) 

320  let  v = w(l)+e*(w(2)-w(l) ) 

325  let  r5  = v*y"(l/7)/sl  ■ • 

330  let  wl  = r5*al 

740  print  " VP.  = "int(*l+.5) 

750  print 

7*0  print  c 

765  ro  to  20  : 

370  rrint  !>0  VALID  VR=  VP  TOO  I, ARC-}.  rDK  FOR 

772  print 

375  ro  to  20  . 

TOO  ren  this  rives  parar.eters  for  the  pwrs 

700  if  hlfi<=P00  then  420  • ' V 

•ioo  Tirint’’HO?  too  larre  for  prof  ran”  , o 

41  f ro  tv.  10  ' ••• 

'■  * 62 
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420  let  h(2)  = h(l)  + 100 
470  lft  e = (hl-h(l)  )/l00 
440  if  h(i)  >=  100  then  "0 
44?  if  v2  > 7.5  then  470 
450  let  kO  = P.206P76 
452  let  V:1  = -9.P66?222e-02 

454  let  k2  = -4. 270571 9e-07 

455  let  k7  = 44. 67751 000 e-05 

459  let  k4  = 0 

460  ro  to  272 

470  let  kO  = P.267247 
47<)  let  kl  = -1.210P524e-01 
4f0  let  k2  = 12.74266e-04 
500  let  k7  = -9.20654P6e-06 
502  let  k4  = 0 
540  50  to  27? 

560  if  h(i)  >=200  then  P60 
556  rerc  COI'.fTANT?  TOR  100  ‘ 
667  if  v2  > 61  then  770 
560if  v2  > 7.5  then  600 
566  let  kO  = P. 29127 
670  let  kl  = -1. 17297 Pe-01 
575  let  k2  = 71 .199 OP e- 06 
500  le  t k7.  = 0 
590  le  t VA  = 0 
696  60  to  £72 
600  let  kO  = B. 29965 
610  let  kl  = - 1 . If 4 777Pe- 0 1 
620  let  k2  = -4.B494f)P5e-04 
626  let  k7  = 55. P7 01e-0 6 
670  let  k4  = -9.16B077Pe-07 
640  r o to  272 

660  if  h( i)  >=  70  0 then  770 
666  reir  CONSTANTS  TV 5 200 
670  if  v2  > 41  then  770 
gap  let  90  = P. 795 227 
690  let  kl  = -1. 471755 6e-0 1 
TOO  let  k2  = 12.744P9e-07 
710  let  k7  = -2.0672771e-07 
720  let  k4  = 16.o7691e-05 
7^0  let  1:5  = -6.R9742P-06 
7^n  let  V.e  = 14.277 14e-Ofi 
750  let  V7  = -1.1 675 0 1 5 e- OP 
760  60  t.o  240 

77p  if  h(i)  >=400  then  PPO 
775  if  v2  > 74  then  570 
TP?  let  kO  = P.41965‘4 
790  let  kl  = -9 .9P27R1 6e-02 
POO  let  k2  = —4. 1B72797 e-07 
PIP  let  k7  = 64.49CP4p— 05 
P2D  let  k^  = -7 ,75R752e-05 
P70  let  k5  = 14.C0969e-07 
P4p  i°  t k.6  = -2.01 70989e-OR 
050  le  t 1:7  = 0 
P70  r o to  240 


ONI  KQ. 


for  ALL  VHS 
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PPO  if  h(l)  >=500  then  PPO 

pp5  wwrrrswKTS  for  ^oo* p kf  fq.  for  all  vks 

. .PPO  if  V2  > 30  then  370  fj 

POO  let  kO  = P.4PP4RP 

PC5  let  kl  = -1.0P6521 le-01  ‘ 

F10  let  k2  = -3. 444574? e-03  ' «. 

P20  let  k3  = ?2.61?06e-05 

PPO  let  k4  = -7.10-!05e-05 

P40  let  k5  = 33.1P013e-07 

P50  let  kfi  = -5.6 6R5 057 e— 0 R 

POO  let  k?  = 0 

PPO  fo  to  240 

PPO  if  h(i)  >=  000  then  1100 

PP2  rem  CONSTANTS  FOR  500* OIF  l.C.for  ALL  VKS 

PP5  if  v2  > 27  then  370 
1C00  let  kO  = P. 525. PRO 

1010  let  kl  = -6.3l20552e-02 

1020  let  k2  = -2.5 622191e-02 

1030  let  1:3  = 54.20447e-04 

1040  letk4  = -5.P2f33Pe-04 
1C-50  let  k5  = 34. R5504e-06 
1000  let  kO  s - 1.02286.46^-0  6 
1070  let  k7  a 11.4432e-0P 
10PC  rc  to  240 

1100  if  h(i)  >=700  then  1210 

1105  reir  CO KSTA  NTS  FOR  500  * ONI.  IC.  for  ALL  VKS 

1107  if  v2  > 25  then  370 
1110  let  kO  = P .5R0222 

1120  letkl  = -1.002711e-0!  -,W 

1130  let  k2  = -P.P17l7.5cft_c3 

1135  let  k3  = 2P.0232P-R4 

1140  let  k4  = -3.502P22.4e-04 

1150  let  1:5  = ?P.02515e-0fi 

1160  let  k6  = -1.0R?6364e-06  . 

1170  let  k7  = 15.4 Iff 7e-0P 
1200  fro  to  240 

1210  if  h(i)  >=P00  then  1320 

1?12  if  v?  > 22  then  370  ■■  . 

1215  rer  001  STAFFS  FOR  70n* OM  J.C.  for  ALL  VN *s 

1220  let  kO  = P.655 9P2 

1225  letkl  = -l.3F7ppfpf._ci 

1230  let  k?  = 14. 262 PI e-03  ‘ 

1235  let  k3  = -4 . 0P2PPP3 e-03 
- 1240  let  k4  = 50. 2 PI  25  e -05  - , , 

1245  let  k5  = -2.57i??3Pe-05 

125-0  let  k6  = 43.'7POC3e-Oc 

1253  let  k7  = 0 t , 

. 1260  eo  to  ?40 

•13.20  ifh(i)  >=500  then  13F0 

1221  if  v2  > 21  then  370  ' . , 

1322  rer.  CONSTANTS  FOR  POO  ■ "=  , 1 ' 

1324  if  v?  > 7.5. then  1350 

,1326  let  kO  = P.6R12P5  * ' \ * • ' e ' 

1330  let  H = -1.1432664^-01  - r» 

, 1335  let  l'_2  = - 1 .78PP6P6e -C3  =‘  ‘ ” 64  « •'  •/  • ° • •• 0 • 
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1340  If  t k.3  = l5.95909e-05 
134?  lft  k4  =0 
1345  po  to  232 
1350  lft  kO  = 12.51342 
1355  let  kl  = -1.516344 
1360  let  k?.  = 17.69944P-02 
1365  let  k3  = -8.900R35e-03 
1370  lft  k4  = 14.00736e-05 
1375  60  to  232 
1380  rens  CONSTANTS  FOR  900' 

13R1  if  v2  > 20  then  370 
1382  if  v2  > 7.5  then  1400 
1385  let  kO  = 8.719654 
1387  let  kl  = -1.215R526e-01 
1390  let  k2  = 12. 03604 e-04 
1392  let  k3  = -1.3R63281e-04 

1394  let  k4  * 0 

1395  go  to  232 

1400  lft  kO  = 13.47289 
1405  let  kl  = -1.971983 
1407  let  k2  = 25 .47267 e-02 

1409  let  k3  « -1. 432.5 115e-02 

1410  let  k4  = 26.40371e-05 
1415  go  to  232 

1430  rem  THIS  SICTION  PKTSF.NTS  PARANF.TIRS  TOR  QWS 

14^0  if  hl<=  900  then  1470 

1450  prinfHOB  is  too  laiy?e  for  propTarc  ’ 

1460  fo  to  10 
1470  lft  h(2)  = h(l)  +100 
1480  le  t e = (hl-h(l)  )/l00 
1490  If  h(i)  >=  100  then  1600 
1494  retn  CONSTANTS  TOR  0 FOB 
1496  if  v2  > 36  then  370 
1600  let  kO  = 8.315159 

1610  let  kl  = -.1060868 

1620  let  k?=  .0005224 
1530  let  k3  = -.000313 

1540  let  k4  = 3.22649e-05 

1550  let  k5  = -1. 23227 e-06 

1660  let  k6  = 1.9f7707p-08 

15W  let  k7  = -1.05880 P-10 

1590  roto240 

1600  If  h(i)  >=200  then  1700 
1604  remCOHSTT NTS  FOR  100'  SPOB 
1606  if  v2  > 35  then  370 
1610  let  kO  = 8.376082. 

1620  let  kl  = -.1042945 

1630  let  k2  = -.0012014 

1640  let  k3  = -3.91l36e-05 

1650  let  k4  = 1.287  57e-05 

1660  let  k5  = -4.9757pe-C7 

1670  let  k6  = 5.77257e-09 

1680  let  k7  = 0 

1695  po  to  240  65  ' 

1700  if  '\)  >=  300  then  1810 
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1704  rem  CONSTANTS  FOR  200*  SHOE 
1705  if  v2>  35  then  370 
1710  le  t kO  = 8.4?0?4 
1720  let  kl  = -1.09473e-01 
1730  let  V2  = 14.622R8e-04 
1740  le  t k3  = -5.9897 92e-04 
17R0  let  k4  = G6.97002e-06 
l^pn  let  k5  = -3.0149459e-06 
1770  let  k6  = 6l.PR22Re-09 
1780  let  Yy  = -4.Rf,6F33e-10 
1C00  ro  to  240 
1810  If  h(i)  >=400  then  1920 
1812  rem  CONSTANTS  KIR  300*  SHOE 
1816  if  v2  > 35  tJien  370 
1820  let  kO  = R. 485315 
1830  let  kl  = -.1031396 
1R40  let  k2  = -.0034114 
1850  let  k3  = .00030R7 
1860  let  k.4  - -1.07267 e-05 
1870  let  k5  = 3.15662e-07 
1880  let  k6  = -5.56646e-09 
1890  let  k7  = 0 
1910  ro  to  240 

1920  if  h(i)  >=500  then  2030 
1922  remCONSTAMS  >DH  400*  SFOE 
192*  if  V2.  > 31  then  370 
1930  let  kO  = 8.576003 
1940  let  kl  = -.1039885 
1950  let  k2  = —.0065768 
I960  let  k3  = .0012382. 

1970  let  k4  = -.0001333 
1980  let  k5  = 8.01387e-06 
1990  let  kfi  = -2 . 34684 e- 07 
2000  let  k7  = 2.5 1295e-09 
2020  ro  to  240 

2030  if  h(i)  >=600  then  2140 

2032  ren  CONSTANTS  FOF.  500'SKOB 

2035  if  V?  > 28  then  370 

2040  let  kO  = 8.643504 

2050  let  kl  = -.1110564 

2060  let  k2  = -.0041901 

2070  let  k3  = .0006644 

2080  let  1:4  = -7. 76848  e-05 

2090  let  k5  = 5.PP6P5e-06 

2100  let  kR  = -2.27C79e-07 

2110  let  K7  = 3.00626e-09 

2130  ro  to  ?4C 

2140  if  h(i)  >=700  then  2250 

2142.  rem  CONSTANTS  TOR  600*  SPOB 

23  4*  IF  V2  > 26  THEN  370 

2150  let  kO  = 8.6856 97 

2160  let  kl  = -.11648.22 

2170  let  k2.  = .1X103634  - 

2180  let  k3  = -.0006169 

2190  let  k4  = R.5754le-05  . , 
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220 0 let  k5  = -4.07263e-O6 
221 0 le  t k6  = 5.66402e-O8 
2220  let  k7  . 0 
2240  ro  to  240 

2250  if  h(i)  >=8C0  then  2360 
2252  rem  COKS  TAKTS  TO R 700*  SHOB 
2255  if  v2  > 25  then  370 
2260  let  kC  = 8. 70744? 

2270  let  kl  = -.1175502 
2280  let  k2  = .0023483 
22P0  let  k3  = -.0O13O54 
2300  let  k4  = .000190? 

2310  let  k5  = -1.15200e-05 
2320  let  k€  = ?.P3079e-07 
2330  let  k7  = -2.447O4e-09 
2350  ra  to  240 

2360  if  h(i)  = 900  then  2470 
2362  rem  CONSTANTS  FOR  800*  SHOP 
23.65  if  v2  > 23  then  370 
2370  let  kO  = 8.73CT28 
2380  let  kl  = -.1151805 
2390  let  k2  = .0021175 
2400  let  k.3  = -.0015218 
2410  let  k4  = .0002654 
2420  let  k5  = -l.P€750e-05 
2430  let  k6  = 6.18015e-07 
2440  let  k7  = -7.20562e-09 
2460  60  to  240 

2470  ren  Constants  for  900  SHOB 
2475  if  v2  > 22  then  370 
2480  let  kO  = 8.793042 
24PO  let  kl  = -.1154885 
25.O0  let  k2  = .0001871 
2510  let  k3  = -.C011O08 
2520  let  k4  = .0002357 
2530  let  k5  = -2.01562.e-C5 
2540  let  k6  = 6.P7520e-O7 
2550  let  k7  = -8.74066e-OP 
2570  60  to  240 
2620  end 
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Computation  of  Probability  of  Damage  to  Point 
Targets  or  Normally  Distributed  Area  Targets 


This  appendix  explains  the  basic  probability  of  damage 
calculation  program  presented  in  Section  VI  for  point  targets 
or  normally  distributed  area  targets  located  either  at  the 
DGZ  or  offset  some  distance  "x."  The  example  program  also 
provides  for  calculating  probability  of  damage  to  a target  when 
CEP  = 0,  i.e.,  when  the  impact  point  is  known.  The  methodology 
outlined  here  may  be  adapted  to  the  calculations  of  probabilities 
of  damage  to  uniform  area  targets  as  explained  in  Section  VI. 


The  probability  of  damaging  a point  target  located  a 


distance  "x"  from  the  DGZ  is 


ITT  *o 

■ / f p<'-: 


where  P^(r) 


Ztt  <rz 


r*t  xz-  Zrx  cos 9 

2<ra  , . 

G rdrdG 


z (r)  = 


. H*/_ 

= ( -~=.  e 
J YztF  c 

- OO 


a = CEP/1 . 1774 


The  probability  of  damaging  a normally  distributed  area 
target  can  also  be  described  by  the  above  equation  if  a is 
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replaced  by  aA  where  ° . 

cepa 

°A  

A 1.1774 

CEPa  = ’''CEP2  + .231*  (P-95)  2 

As  expressed  above  P = f (3 ,WR,x,o ,P-95) . In  computing 
P.D.  slightly  different  parameters  are  input.  They  are 
3=0^=  damage  sigma  (.1,  .2,  .3,  .5) 

where  3 = sqr  (-ln{l-s2)) 
wl  = WR 

xl  = x “ offset  distance,  which  may  be  zero 

c = CEP  = measure  of  delivery  error.  (When  P.D.  from  GZ 
is  desired,  input  c = 0) 

o = CEP/1 . 1774 

ca  = CEPA  = adjusted  CEP , 

cepa 

and  aA  = — - where  a l describes  the  variance  of  the 

A 1.1774  A 

joint  distribution  of  delivery  error  and  target 
density. 

tl  = P-95  = radius  of  circle  which  encompasses  95%  of  the 
target  being  considered. 

In  the  example  program  CEPA  = ca  = *1 c2+ .231* (tl) 2 . 

Except  for  the  case  of  calculating  P^fr)  from  a GZ  to  a point 
target,  all  variables  are  standarized  by  division  by  aA  = 


1.1774. 


If  the  target  is  a point  target  then  tl  = 0 and  ca  = c. 
If  the  probability  of  damage  to  a normally  distributed  area 
target  from  a GZ  is  desired  then  c = 0 and  ca  = 231* (tl) 2 . 

A CEP  of  zero  denotes  a situation  in  which  the  impact 
point  is  known  and  it  is  necessary  to  compute  the  probability 
of  damage  to  a target  located  "xl"  distance  from  the  impact 
point.  This  probability  is  calculated  by  letting 
z(xl)  = 1/8  * In  (wl*exp (-3 2 ) /xl) 


The  erf  may  be  approximated  as  shown  in  Section  VI. A. 3. 
The  approximation  is  fna(v)  in  the  example  program.  If  z>3.87 
then  P(j=l  and  if  z < -3.87  then  = 0 since: 


and 


< .0001 


71 


If  Jzj  < 5 x 10“ 6 then  P<j  = .5.  This  is  to  avoid  an 
exponential  underflow  on  the  computer,  since 


u 


erf  (vlD  “ ^na  (v)  ~ 


I 


where  V = 


lal 

nr 


and  D = I + e,v  + ezvz+-  e3v3+  ^v^ejV’te^v* 

When  J z{  < 5 x 10"6,  D18  will  contain  negative  exponents  too 
large  for  the  computer  on  which  these  example  programs  were 
written  to  handle  and  an  exponent  underflow  message  would  be 
printed.  To  allieviate  this  proM®:;,  is  set  equal  to  .5 
when  I zf  < 5 x 10"6. 

When  deliveiy  error  and/or  target  density  must  be  con- 
sidered, all  parameters  are  standardized  by  division  with 
°A  “ CEPa/1.1774 

w3  = — = 1.1774  * wl/ca 
°A  a 

'«  x = ~ = 1.1774  * wl/ca 


r = (distance  from  impact  point  to  target) /oA 
The  probability  of  damage  P is 


P = ^ -f(r)  dr 

‘ -£(r%-xz) 

where  -f(t-)  — f^Cr)  r tS 


2ir 


j-  ( 


r*  cos  © 

e de 
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Since  it  is  necessary  to  evaluate  f(r)  only  as  long  as 
f (r)  is  significantly  greater  than  0, 


cO  D 

^ -f(r)  r = ^f(r)di 


where  f (r)  z 0 whenever  r<a  or  r>b  . 

Wherever  r = 9.2  5 * WR  * damage  sigma  (cr^  = .2,  .3,  .4, 
or  .5),  Pd(r)  £ .0005,  and  it  is  unnecessary  to  evaluate  the 
integral  beyond  the  point  b = 9.25  * WR  * . 

Impact  points  are  assumed  to  be  normally  distributed 
about  the  DGZ . The  probability  that  a weapon  would  impact 
farther  than  4a  from  the  DGZ  is  given  by 


zn 


oo  _*/ 

r -rA  i 

f re  dr  a oooc5 


Since  the  DGZ  is  located  "x"  distance  from  the  target,  there 
is  no  significant  contribution  to  P.D.  when  r < (x-4)  or 
r > (x+4)  . Therefore,  a = max  (0,x-4)  and  b = min  ( 10  *WR*a(j  ,-x+4 ) 

The  integral  \ ■f^’r)dr  may  be  evaluated  using  a 10  point 
Gauss -Legend re  quadrature  formula  i.e., 

b 

io 


where 


(f(r)dr  = E 


rt  = . 5[(b-a)  + b + a]  . 
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The  base  points,  , are  symmetrically  placed  with  respect  to 
the  midpoint  of  the  interval  of  integration.  The  weights,  w-^  r 
are  the  same  for  each  symmetric  pair  of  Zj_ . It  is  customary 
to  designate  for  the  pair  of  symmetric  points  nearest  the 
midpoint,  z ^ for  the  next  symmetric  pair,  etc.  The  notation 
z±n  = J,  w±n  = K means  zn  = J , z_n  = - J , wR  = K,  w_n  = K. 

10-Point  Formula 


±1 

a 

0.14887 

4339 

W±1 

.29552 

42247 

±2 

= 

0.43339 

53941 

w±2 

= 

.26926 

67193 

±3 

a 

0.67940 

95683 

W±3 

= 

.21908 

63625 

±4 

a 

0.86506 

33667 

W±4 

a 

.14945 

13492 

±5 

= 

0.97390 

65285 

w±5 

a 

.06667 

13433 

At  each  r,  z = 1/0  ln(wl*exp(-02)/r)  and  pi  = P<a(r)  is 
calculated  as  previously  explained  using  fna(v).  If  h=0 , 
i.e.,  x=0,  then  f (r)  = pi*r*exp 1-rYz) 

If  h = (x*r)  < 3.75  then  j = (h/3.75)2 

and  f (r)  = pl*r*exp(-(x2+r2)/2) *fnb ( j ) 

I 

where  f nb  ( j ) . j exp(xrcose)  d0 

o 

If  h = x*r  > 3.75  then  j = 3.75/h 

and  f(r)  = pl*r*exp  {- (x-r)^L)  * fnc  (3  ) 

J & 


where  fnc (3 ) *exp (xr) 

& 


= j — exp (xrcosO) d9 


f (r)  is  then  weighted  by  wj_  and  summed  into  f according  to 
the  Gauss-Legendre  quadrature  formula 


P = . 5 (b-a)  *f 
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SYS  TIP  ? . ■ ■ *>■ 

basic 

old  or  new-oiri  newprob  c 

FIIJ  nevprob  CLASSIFIED  OUT 

ready 

♦list 

10  re  it  CALCULATE  P.D.  TO  POIKT  TARGETS  10-PT.  Gaussian-Legendre 

20  dim  7(5)  ,v(5) ,f (5,2 ) 

40  data  .0^6*2308, ’.0422  R2.01,  .00927053,  .00015201 

45  data  .0002  7657  , .00004306 

£0  def  fna(v)  = l-l/(l+v*(el+v*(e2+v*(e3+v*(e4+v*(e5+v*e6))))))  16 
60  read  bl ,h2,b3,b4 ,h5,b6 
70  data  3.5156229,  3.0899424,  1.2067492 
80  data  .2659732,  .036076R,  .0045813 
SO  def  fnh(v)  =l+v*(M+v*(b2+v*(b3+v*(b4+v*(b5+v*b6))))) 

100  read  c0,cl,c2,c3,c4,c5,c6,c7,cfi 

110  data  .398942 2P, .01 32P592,  .00225319,  .00157565 

120  data  .00916281,  .02057706,  .02635537 

130  data  .01647633,  .0039^377 

140  def  frc(v)  = c0+v*(c1+v*(c2-t#(c3-7#(c4-v*(c5-v*(c6-'s^c7-v#c8)  )))))) 

150  for  I = 1 to  5 
160  read  z(i),v(i) 

170  next  i 

180  data  .148874339, .295524225  ,.433395394,  .269266719 

185  data  .679409568,  .2190863©,  .865063367,  .149451349 

190  data  .973906529,  .066671344 

200  nrint"lNPUT  SIGMA-  (decimal  form) " 

230  input  s 

220  print"lNPUT  VP,  X,  CEP,  P-95” 

230  input  vl,xl,c,tl 

240  let  h = sqr(-lof (1-s*'?.)) 

245  let  c = sqr(c"2+.231*tl  2) 

250  If  c > 0 then  430  < 

260  if  xl  >0  then  290 
270  let  p = 1 
280  60  to  R60 

290  let  z = (l/b)*lof((vl*exp(-b"2))/xl) 

300  if  z > 3.87  then  270 

310  if  abs(z)  < 5e-06  then  390 

320  if  z < (©.87  ) then  *10  * • 

330  let  u = abs(z)/sqr(2)  « 

340  let  e = fna(u) 

350  let  t»  = .5  + .5*abs(z)/z*e 

360  60  to  860  * 

370  let  p = 1 
380  fo  to  860 

390  let  p = .5  ■ . 

400  60  to  860 

410  le  t p = 0 76 

4?0  60  to  86.0 

S 27050 
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470  1ft  = 1.1974*wl/r 
440  1ft  X = U774*xl/c 
450  If  t f = 0 
460  If  t bP  = P.25*w7*s 
47 0 If  bP  < 4 thfn  WO 
460  1ft  a?  = x-4 
4PC  if  aP  > 0 thfn  510 
500  If t aT  = 0 
510  let  bP  = x+4 
520  eo  to  540 
WO  let  aP  = 0 
540  for  i = 1 to  5 
550  for  n s 1 to  2 

560  1ft  r = .5*(z(i)*(-l)  n*(bp-af-)  + bP+aP) 

5”P  If t z = (l/b)#lof ((w5*exp(-b"2) )/r) 

500  if  z >5.57  thfn  650 
5P0  if  abs(z)  < 5e-06  thfn  670 
600  if  z < (-7. 67)  thfn  6P0 
61 C 1ft  u = abs(z)/sqr(? ) 

620  1ft  f = fna(u) 

WO  1ft  pi  = .5  + .5*abs(z)/z*e 

64 0 ro  to  710 

WO  1ft  pi  = 1 

660  fo  to  7 10 

670  1ft  pi  = .5 

6«0  “ 0 to  710 

(TO  If  t pi  = 0 

TOO  fo  to  7i0 

710  lft  h = x*r 

7?0  if  h = 0 thfn  600 

770  if  h > 7.75  thfn  770 

74  0 lft  J = (h*h)/( 7.75*7. 75)^ 

750  lft  f,A,n)  = plvr*exo(-(x~2+r'2)/2.)*fnb(J) 

7 fp  ro  to  620 
770  let  .1  = 7.76/h 

7P0  lft  f(i,r)  = pl*r*exp(-(x-r)’2/2)*fnr(  j)/sqr(h) 
TfO  ro  to  P20 

POO  lft  ffi.n)  = pl*r*f xp(-r"2/2) 

71  0 p n to  620 

P20  lft  f = f*v(i)*f (i,r.) 

670  nfxt  r. 

P40  nfxt  i 

.TO  lft  p = ,5*(hP-aP)*f 

P6 o "rir.t  ir.ti  1000*r+.5)/1000 

77  o print 

«70  rv  to  270 

<Tr  *“nd 

ppipv 
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7 

♦run 


HFUT  SIGI-A-  (decimal  form) 

•> 


IK?t’T  VP,  >,  CF.P,  F-9S 
7 

3500, 2500, 1000,0 
.CPC 


•? 


2250,0,1500,500 

.CP? 


.WOO, 1000 ,500,0 

,?P9 


1500,0, ,f 

.*?C 


7 

1500,0^00,0 

.POP 

7 

if oo  for, looo 
.per 


•> 
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Computation  of  the  Offset  Distance  for 
a Specified  Probability  of  Damage 

Situations  arise  when  it  is  desirable  to  know  the  maximum 
offset  distance,  x^  at  which  an  aim  point  may  be  placed  and 
still  maintain  a specified  probability  of  damaging  a target. 

It  is  possible  to  find  this  maximum  offset  distance  by  an 
iterative  process  of  increasing  the  offset  distance  and  cal- 
culating P.D.  until  the  lowest  acceptable  P.D.,  Pm,  is  found. 
However,  this  method  is  slow  and  would  be  much  too  time  con- 
suming if  many  cases  had  to  be  considered. 

Equations  have  been  developed  which  yield  a rapid  calcula- 
tion of  xM  when  45%  < < 95%.  Two  equations  are  given  for 

each  damage  sigma,  an  equation  of  a curve  is  used  for  smaller 
values  of  WR/CEP  and  a linear  equation  for  larger  values. 
Constants  in  the  equations  are  functions  of  Pm. 

Let  x’j  demote  the  offset  distance  such  that  P (x^)  = Pm 
and  let  xM  be  the  offset  distance  as  determined  by  the  equations. 
The  following  statements  of  accuracy  generally  hold. 

Ixm_xmI  ^ 0.10  CEP 

In  the  range  of  the  equations  of  a curve,  i.e.  for  smaller 
WR/CEP  : 

|Pm-P(xM}|  S .02 

Appendix  F 
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In  the  range  of  the  linear  equations,  i.e.,  for  larger  WR/CEP: 
|pnri><>‘M>i  S -01 

There  are  two  exceptions  to  this  general  accuracy,  both 
occur  at  od  = .50.  When  Pm  is  in  the  range  of  60%  then 
max  Jpm-p (xM)|  < 2.5%.  When  P is  the  range  of  50%  then 
max  |Pm"p(xM)|  <5%.  If  greater  accuracy  is  required,  it  is 
suggested  that  the  equations  be  used  to  derive  xM,  then  xM 
used  as  an  initial  starting  value  in  the  P.D.  computer  program, 
increasing  or  decreasing  xm  as  necessary  to  find  x^» 

The  following  equations  are  used  for  determining  xM  as 
f(od,  WR,  CEP,  Pm) . For  SIGMA-20  targets,  when  WR/CEP  < 3,  then 

xM  = CEP*  A*  (wl-wO) B 

where  wl  = WR/CEP 

wO*  = .58  + 1.7*Pm2 

A = 5.102  - 14.9435*Pm  + 

22.515‘Pu,2  - 11.674*Pm3 
B = -5.218  + 27.37*Pm  + 

-41.766*Pm2  + 20.664*Pm3 


*w0  indicates  the  smallest  WR/CEP  which  can  achieve  Pm  at  the 
DGZ . Therefore  in  the  sample  program  if  wl  < wO  a statement 
is  printed  saying  "P  cannot  be  achieved  with  this  weapon." 
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When  WR/CEP  > 3 


where 


For 


waere 


When 


wnere 


xM  = CEP* ( A + B*wl  ) 
wl  = WR/CEP 

A = -(.0153  + 1.184*Pm4) 
B = 1.195  - . 456*Pm 

Sigma-30  targets,  when  WR/CEP  < 3.5, 
xM  = CEP*(a+B*  (wl-w0)+C*(w-w0)  2) 
wl  = WR/CEP 

wO  = (1.522  - 1.2058*Pra}-1 
A = 1.146  - 1.85 *Pm+ .975 *Pm z 
B = 1.4 

C = -.414  + 1.126*Pm  - 1.14 *Pm 2 
WR/CEP  > 3.5, 

= CEP* (A  + B + wl) , 
wl  = WR/CEP 

A =-(.00662  + .93  * Pm4) 

B = 1.256  - -678*Pm 
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For  SIGMA-40  targets,  when  WR/CEP  < 3.5, 

xM  = CEP*(A+B(wl-wO)  + C(wl-wO)2) 

where  wl  = WR/CEP 

wO  = (1.467  - 1.217*Pm)_1 
A = Pm/ (-2. 73  + 8.033  *Pm) 

B = 1.7  - .4115  *Pro 
C = -{.508  - 1.184*Pm  + 1.118*Pm2) 

When  WR/CEP  >3.5, 

xM  = CEP  * (A  + B * wl) 

where  wl  = WR/CEP 

A * .026  - l.Ol*?^ 

B = 1.234  - .785 *P 

m 

For  Sigma-50  targets,  when  WR/CEP  < 5, 

; ‘ . «1  ,,  r rr 

« XM  = CEP*^A+B (wl-wO) +C* (wl-wO ) 

where  wl  - WR/CEP 

wO  = (1.34  - 1.15*Pm)_1 
A = .28 

B = 2.152  - 1.524*Pm 
c C = -(.0252  + . 35*Pm  - . 33*Pm2 ) 
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When  WR/CEP  > 5 


xM  = CEP* (A  + B * wl) 

where  wl  = WR/CEP 

A = .113  - .864‘Pjn4 
B = 1.17  - .858  * Pm 
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3 

♦list 

03  rem  CALCULATES  MAX.  OFFSET  FOR  MIN.  REQUIRED  P.D. 

05  print "INPUT  S” 

10  Input,  $ 

15  print”lNPUT  P,  WR,  CEP" 

20  Input  p.w.cO 

25  let  wl  s w/cO 

30  If  s b .2  then  200 

35  if  s = .3  then  300 

40  If  s * .4  then  400 

45  if  s = .5  then  500 

50  print"nonappropriat,e  sistn&" 

55  *o  to  10 

200  if  vl>  3 then  250 

210  let  WO  * ,5ft+1.7*p*> 

212  if  wl<  wO  then  500 

215  let  a = 5. 102-14, 9435*D+22.515*p~2-lt.674*p'*3 
220  let  b = -5.210  + 2?.37*p  -41.766V*?.  +20.664*p~3 
230  let  xl  = a»(wl-wO)"b 
240  *0  to  600 

250  let  a * -(.0153  * 1.1R4W) 

255  let  b a 1.195  - ,456*p 
260  let  xl  = a+h*wl 
2.65  fro  to  600 
300  if  wl  > 3.5  then  350 
305  let  vO  = (1.522-1 ,205fl*p) ~(-l) 

307  if  wl<  wO  then  500 

310  let  a a 1 . 146-1 . 05* D+ . 975 »p~2 

315  let  b s 1.4 

320  let  c = -.414  + 1.126*p  -1.14V*2 
025  let  xl  = a+b»(wl-wO)  + c*(wl-w0)"2 
330  *o  to  600 

350  let  a a -(.00662  +.93*p~4) 

355  let  b = 1.2.56  -,67R*p 
360  let  xl  a a+b*vl 
365  r,o  to  600 
100  if  wl  > 3.5  then  450 
405  let  wO  = ( 1.467-1. 2 17*p)'*(-l) 

407  if  wl  < wO  then  500 

410  let  a = p/(-2.73  + fl.033*p) 

415  let  b = 1.7-.4H5*p 

420  let  C = -(.506  -l,164*p  + l.llR*p~2) 
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425  let  xl  = R+h*(wl-wn)  + c*(vl-vOV'2 
430  ko  to  «in 
4.00  let  a = .026  -1.0l*p~4 
455  let  h = l.?34-.78T.*p 
460  let.  >l  = &4-h*wl 
4P0  PO  to  600 
500  if  wl  > 5 then  .550 
.605  let  wO  = ( 1 .34—1 . 15*o) '(— 1 ) 

507  if  wl  < vO  then  500 

510  let  a = .20 

51.'  let  b = 2. 152-1. 524 *p 

520  let  C - -(.0252  +.35*p  -.33*p'%2) 

5£5  let,  xl  » a+b*(wl-wO)  +c*(wl-v0)"2 
530  ao  to  600 
550  let.  a = .113-. 064 
55.)  let  b = 1.17-.05fl»p 
560  let  xl  = a+t*wl 
565  »o  to  600 

500  printp.w,  p cannot  be  achieved  with  this  weapon" 

505  e o t.o  20 

600  if  xl5  0 then  fin-, 

603  let,  xl  = 0 
605  let,  x = xl*cO 
610  print"  "x 

620  print. 

630  po  to  20 
090  end 

r^adv 

•run 


IH°UT  S 

? .2 

T),  WP,  C2P 
7.8,4500, 1000 

3235.634 

7.9,5500,1500 

3127.116 

7.75,10000,2500 

7555.187 

7.6,10000,2000 

8« 76 .507 

? .5,2000,1000 

1794. ?R9 
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?*.95,?oho,iono 

.95  ?^nn 

p cannot  be  achieved  with  this  weapon 
?. 05 ,2000,1000 

503.170? 


•run 


INPUT  S 

?.3 

INPUT  ?,  KP,  cr.5 

? .75,4500,3000 

.7 5 450(0 

p cannot  be  achieved  with  this  weapon 
?. 75 ,9000, 3000 

5515.il? 

? .05,5000,1500 

?457. 107 
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Computation  of  Probability  of  Damage  to  Circular 
Area  Targets  with  Uniform  Density 


To  compute  probability  of  damage  to  circular  area  targets 
with  uniformly  distributed  elements,  concentric  circles  are 
drawn  about  the  DGZ  (or  GZ  if  CEP  is  0)  with  radii  rj_  = r^-i  + 
Ar  where  Ar  is  a small  increment.  The  i^*1  annulus  is  defined 
as  lnrx-l  it  r*,  • The  probability  of  damage  to  any  target 

element  located  in  the  ith  annulus  is  Pj_.  If  ai  is  the  area 
of  intersection  of  the  ifc^  annulus  with  the  target  then  the 
average  probability  of  damage  to  the  uniform  circular  target  is 


Pi.  au 

P=.-L 

Sa- 

i=I 

where  TR  = target  radius  f 


E«,  = tttr' 


and  I is  such  that  a circle  of  radius  rj  about 
the  DGZ  encompasses  the  target. 


WR,  target  radius,  offset  distance  cf  the  center  of  the 
tarqet  to  the  DGZ  and  CEP  are  input  into  the  example  program 
as  wl , tl,  xl , and  c.  Computation,  however,  is  based  on  the 
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parameters  w3,  t3,  and  x3.  If  CEP  = 0 then 


w3  = WR  = VI 


If  CEP  > 0,  then 


t3  = tl  = target  radius 
x3  = xl  * offset  distance 


w3  = 


1 . 1774  *wl  WR 


t3  = 


1.1774*tl  TR 


x3  ~ - . 1774 *xl  _ offset  distance 
c a 


The  radius  of  the  itla  concentric  circle  about  the  DGZ  is 
r(i)  = rO  + i*kl.  rO  is  the  initial  value  of  r(i)  and  is  0 
if  the  DGZ  is  located  inside  the  target  area.  Referring  to 
figure  G-l , if  the  DGZ  is  located  outside  the  target  area 
then  rO  = x3-t3.  In  the  example  program  kl  = Ar  = t3/10 . 


ro-o 

DGZ  Inside  Target 


Figure  G-l 


ro=x3- t3 

DGZ  Outside  Target 


%■ 


If  the  intersection  of  the  i^h  concentric  circle  about 
the  DGZ  with  the  target  is  s(i),  then  the  intersection  of  the 
ith  annulus  with  the  target  is  a(i)  = s(i)-s(i-l).  (See  figure 
G-2) 


Figure  G-2 


The  program  uses  the  notation  x = r(i)  -.5«kl,  where  x is  the 
average  offset  distance  of  target  elements  in  the  ith  annulus. 

The  average  probability  of  damage  p(i)  to  target  elements 
in  the  ifc^  annulus  is  computed  as  follows. 

If  the  CEP  = 0,  then  p(i)  = Pd  (x)  = 1 C di j 

where  z (x)  = -jjp-  Jim  — - ) 

If  CEP  > 0,  then  p(i)  = ^ ^ r exp  (- . 5 (r  2+x2-2rxcos0  ))drd£ 


If  CEP  > 0,  then  p(i 


exp  ( - . 5 (r  2+x-2rxcos0  ))drd0 


The  probability  p(i)  is  multiplied  by  a(i)  and  this 
product  is  summed  for  all  i. 
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This  process  is  continued  by  incrementing  r(i)  , computing 
s(i),  a(i),  and  p(i)  for  each  r(i),  and  summing  p(i)*a(i)  until 
r(i)  = t3  + x3,  i.e.,  until  the  concentric  circles  about  the 
DGZ  cover  the  total  target  area. 

The  average  P.D.  to  the  target,  P,  is  then 

S p a)aa) 



^ ( t3)a 

The  area  of  intersection  of  two  circles,  s(i),is  cal- 
culated as  follows.  If  the  DGZ  is  inside  the  target,  then 
r (i)  < t3  = x3. 

In  this  case  s (i)  = tr  r{i)2  (See  figure  G-3) 


Figure  G-3 


For  a DGZ  inside  or  outside  the  target  where  r(i)  is  greater 
than  (t3  - x3)  (See  figure  G-4) , 
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the  following  notation  is  used 

dl  = the  cartesian  coordinate  in  the  horizontal  direction  of 
the  intersection  of  the  ith  circle  about  the  DGZ  with 
the  target,  based  on  a coordinate  system  centered  at  the 
DGZ. 

_ rar*- t3*  + X3* 

Z « *3 

d2  = the  cartesian  coordinate  in  the  horizontal  direction  of 
the  intersection  of  the  i^-*1  circle  about  the  DGZ  with 
the  target,  based  on  a coordinate  system  centered  at  the 
target  center. 

= dl-X3 

yl  = y2  = the  cartesian  coordinate  in  the  vertical  direction 
of  the  intersection  of  the  ifch  circle  about  the 
DGZ  with  the  target.  It  is  the  same  if  the  coordinate 
system  is  centered  ac  the  DGZ  or  at  the  target  center. 

yl  = ->/ra)a-dli 

yZ=  Vt3a  -d 

al  = atn  (yl/dl)  "J 

f where  atn  is  the  arctangent  function 
a2  = atn  (y2/d2)  J 
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In  the  example  program,  if  yl/dl  or  y2/d2  is  negative  the 
computer  places  the  arctangent  in  the  fourth  quadran:.  Since 
an  angle  in  the  second  quadrant  is  desired,  n radians  =ire 
added.  Therefore, 

if  dl  < 0 then  al  = atn  (yl/dl)  + it 

if  d2  < 0 then  a2  = atn  (y2/d2)  + u 

if  dl  or  d2  = 0 then  al  or  a2  is  tt/2  radians 
s (i)  = si  + s2 

where  si  = r (i) 2* (al-. 5sin (2 (al) ) 

s2  = t32*ii  - t3 2 * (a2-. 5sin  (2*a2) ) 
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RFAPY 

♦list 

Of;  dim  x3  (26)  ,r(25).a(26)  ,s(26)  ,f(6,2) 

06  dim  z(6),v(6) 

1C  read  el,e2,e3,e4,e6,e6 

IF  data  .070623  0784,  <■04  228201?-^  . .009270672 

20  data  .00016?  0143,  .0002766  67?.,  .00004  3 0638 

2?  def  fna(v)  = l-l/(l+v*(el+v*(p2+v*(e3+v*(p4+v*(e5+v*e6)}))))  16 

30  read  hi ,h?,fc3,b4 ,h6,h6 

3P  data  3.6166229,  7.0899424,  1.20674P2 

4C  data  .2669732,  .0360768,  .0046613 

46  daf  fr.b(v)  = l+v*(tl+v*(h2+v*(h3v/*(te+v*(b5+v*te) }) )) 

60  read  cO ,cl,o2,c3,<H,c5,c6,c7,e8 
66  data  .39094228  , .0132P6P2  , .00226319,  .00167666 
60  data  .00916261,  .02067706,  .02636637 
data  .01 64 '76333,  .00392377 

TO  def  fnc(v)  =cO+v*(cl+v*(c2-v*(c3-v*(c4-vMc.Vv*(c6~\r*(</?-v*cPJ)  ))))) 
■>1  for  12  = 1 to  6 

72  read  z(i2),w(i?.) 

73  next  i2 

74  data  .148874339,  .296624326  , .4333P6394,  .269266719 

76  data  .679409668,  .219066363,  .8650633f?',  .149461249 

76  data  .PP390662P,  .066671244 

7p  nrint"  irput  SIGMA  - (decimal  forn)" 

PO  Input  s 

66  print"  input  VR,TF,X,CFP" 

PO  input  vl,tl,xl,e 

96  let  l ~ sor(  -lor(l-s"? ) ) 

100  If  c =0  then  126 

106  let  x2  = 1.1774*xl/r 

110  let  v3  = 1.1774*wl/c 

116  let  t3  = I.1774*tl/C 

120  eo  to  140 

126  let  x3  = xl 

13C:  let  w3  = vl 

136  let  t3  = tl 

140  let  r2  = 0 

146  le  t i = 0 

160  le  t 11  = t7/J0 

166  if  x3  > 0 tier.  170 

160  rn  sub  320 

166  ao  to  90 

170  if  x3  % t.3  then  190 

176  ren  PC 7 11,610}.  TAKf-FT 

1°0  let  rO  = 0 

1°6  ro  to  ?Of 

190  rem  IC  7 Pin  61 ’'I  VAT-KT 
196  let  rO  = r~-t7 

200  rer  rf  i)  IS  Tl:?  KAFKiS  ft?  7:-;T  l*7F  true  ?n  KIP  WHUi.  A.tVi  r.u  iu/. 

206  let  1 = i+J 

210  let  rfil  = rO  +i*):l 

216  if  x3  +r(i)  >t3  then  236  93 
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220  ren  ANNULUS  CONTAin>!1>  INSI1)>:  TAVC-FT 

225  1ft  s(i)  =3.14l5p*r( i)*‘2 

270  eo  to  255  ‘ 

275  ren  ANNULUS  KOT  Cfl.TAIKFD  INSIDK  TARGH* 

240  if  r(i)  < t3+x3  then  250 
245  1ft  r(i)  = t3*-x3 
250  po  sub  685 
255  lft.  a (i)  = s(i)  - s(i-l) 

250  rfn  x IS  TON  AVI  WTO.  OFFSrT  I'lsmCi  C>  7AKCKT  HIM  MS 

261  P1M  I 

255  1ft.  x = r(l)  -.5*1:1 
270  50  Slit  385 
275  lft  p?  = p2+p*ft(i) 

2°C  if  p.ts(d l-r(i))<  .f]  tier.  2f0 
285  ro  to  205 

290  lft  P2  = p2/(7.l4l59*t5'2) 

700  print  usirr  705  ,o2 
705: 

710  print 
315  r o to  90 

320  rfn  DC  7.  AT  CTO  IP  O?  TARCJ.T 

325  for  i = 1 to  10 

330  lft  r(i)  = i*kl 

335  lft  x = r(i)  -,5*kl 

740  lft  s ( i ) = 3.14i59*r(i)*2 

345  let  a(i)  = s(i)-s(i-l) 

350  po  sub  385 

355  let  p2  = p2  +p*a(i) 

35.0  next  i 

355  let  p2  = p2/(3. 14159*t3  2) 

370  print  usinp  305, p2 
3R0  return 

365 rem  CALCULATFS  PD 
3P0  if  c = 0 then  510 
405  let  f = 0 
410  let  1)9  = 9.25*w7*s 

415  if  b9<  4 then  433 

420  let  aP  = x-4 

423  if  a9  > 0 then  428 
425  If  t aS  = 0 

428  let  If  = x+4 


4.70  r o to  435 
473  le  t aP  = 0 
435  for  i2  = 1 tc  5 
476  for  r,  = 1 to  2 


440  let  r = ,5*(  7,(i2)*(-l)  r.*(hr— aP)+hP+aP) 
445  let  z = (l/b)*lop ( (w3*exp(-b"2)  )/r) 

45.0  if  z >7.8.7  then  4 85 
455  if  abs(z)  < 8e-06  then  495 
480  if  z < (-3.87)  then  505 
485  let  ii  = abs(z)/:ar{?) 

4?0  let  e = fna(u) 

475  let  pi  = .f:+.6*ai.s(z)/z*e 

40P,  /.-p  to  51f, 

485  le  t Pi  = 1 94 
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490  fo  to  515 
40S  let  rl  = .5 
500  ro  to  515 
505  let  pi  =0 
510  ro  to  515 
515  let  h = x*r 
5 20  if  h = 0 then  570 
522  if  h > 5.75  then  535 
525  let  J = ( h*h )/ (3 .75*3.75 ) 

527  let  f(i2,n)  = pl*r*exp(-(x  2+r  2)/2)*fnh( j) 

530  ro  to  590 
!Rf  let  J = 3.75/h 

540  let  f(12,n)  s pl*r*exr-(-  (s  r)  '2/2)*fnc( J)/sqr(h) 
545  ro  to  575 
550  ro  to  590 

570  let  f(i2,n)  - pi*r*exp(-r”2/2) 


590  let  f = f+w(12)*f (12,n) 
592  next  n 
•95  next  i2 

POO  let  p = ,5*(b9-a9)*f 


605  r o tc  6P0 

QO  ren  CALCULATES  Pd  FP.0V  GZ 

615  let  z = l/h*lor((vl*exp(-b"2) )/x) 

C20  If  z > 3.97  then  655 
625  if  abs(z)  < 6e-06  then  665 
650  if  z < (-3.R7  ) then  675 
6t5  let  u s abs(z)/s<ir(2) 

6*0  let  e = fna(u) 

645  let  p = ,5+.5*abs(z)/z*e 

05  0 ro  to  680 

655  le  t r = 1 

660  ro  to  69,0 

665,  let  p = .5 

670  ro  to  680 

675  le  t p = 0 

690  re  tn  rr. 

695  ren  CALCULATIS  Il-T} RSIOTION  V11N  OKI  CIRCU  IS  NOT  COmiKLD 

690  retr  IhSlPr.  TPh  OTHER 

695  let  dl  = ( r ( i )’  2 -ttf  2 +x3  2)/(2*x3) 

70 0 let  yl  = sqr(r( i)"2  -dl~2) 

T5  let  d?  = dl-x3 

710  let  y2  = sqrf t3  ’?-d2  *2) 

715  if  d2  =>  0 then  730 

720  let  a2  = atn(v2/d2)  + (355/113) 

725  ro  to  750 

730  i?  d2  > 0 then  745  , = 

■•35  let  a2  = (3.55/113)*. 5 


^40  ro  to  750 
745  let  a2  = atn(y2/d2l 
750  if  dl->  0 then  7g5 
755  let  al  = atn(yl/dl ) + 
760  m to  7a5 


(3,65/113) 
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765  if  dl>  0 then  780 

770  let  al  = (355/ll3)*.5 

775  60  to  7R8 

780  let  al  = atn(yl/dl) 

785  let  si  = r(i)"2*(al-.5*sin (2*al) ) 

700  let  h = a 1*1 80/3. 14159 
^P5  let  hi  = a2*l 80/3. '1159 

800  let  s2  = t3"2* (355/113)- t3''2*(a2-.5*sin(2*a2)) 

805  let  s(i)  = sl+s2 

810  let  a(i)  = s(i)-s(l-l) 

815  return 

825  end 
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input  SJGKA  - (decimal  form) 

? 

.2 

input  VH.TP.X.CPP 


3500, 1500,4000, 1R00 


5000,2000,4000,1500 


.2pa 


.615 


1500,450,500,2000 


,299 


/ 


1500,45(10,0,1000 
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DISTRIBUTION  LIST 


DOD  AND  JOINT  ACTIVITIES 

A01C 

l 

AFS  COL 

a:i5 

2 

OARPA 

A02A 

A 

CEF  NUCLEAR  AGCY 

AOA5 

2 

DC  A 

AOA8 

2 

OEF  CIV  PREP  AGCY 

A^ 

l 

NAT  WAR  COLLEGE 

A995 

13 

NKCSSC 

AC96 

3 

NEACP 

AC57 

2 

CIA  GRP  ANMCC 

A 1 1 5 

l 

CASO  1 ISA) 

A 125 

2 

CSO  (DD  R C El 

A 17* 

6 

WSEG 
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l 

DMATC 

A219 

l 

OMAHC  CODE  MP2 

A220 

300 

DMAHC  DEPOT  PHIL 

A221 

15 

DMAHC  DEPOT  CFLD 

A 300 

8 

JCS 

A 301 

1 

JCS/J-5  STRATNEGBR 

A 330 

l 
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A353 

9 
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A 51 2 

l 
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1 
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1 
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l 
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B61 2 

l 
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B65 1 

l 
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2 
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1 
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25 
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B825 

1 
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U.S. 
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1 
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' 
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1 
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( 
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* 
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C523 
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C 53  8 
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WHITE  SANDS  MSL  RG 

C 535 

1 
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t 
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i 

C 59 1 

2 
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1 
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1 

INFANTRY  CTR 

j ■ 
1 

C6C5 

2 
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1 

C61  5 

6 

MIA  REDSTONE 

i 

C 632 

2 
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i 

C 635 

l 
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; 

C 683 

l 
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C 763 
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1 
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U.S. 

NAVY 

i 

oone 

2 
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A 

D033 
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? 

D"42 

l 

NAVFACENGCOM 

0151 
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0153 

l 

PAC  MSL  RANGE 

0155 

1 

NAVAIRDEVCEN 

0202 

1 

NAV  WAR  COL 

0217 

1 

NA VWPNSCEN 

’* 

0 22C 

l 

CNR 

0244 

1 

NAVOROSYSCOM 

0246 

l 

NAVWtPSLAS  DAHLGRN 

; 

0247 

2 

NOL  WHITE  OAK 

0 24  p 

4 

NAVSHIPSYSCOM 

0256 

1 
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0369 

l 
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* 

D496 

1 

OIR  SSPO 

0561 

l 
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D947 

1 

NAVINTCOM  34 

0058 

1 

OP  603 

D28C 

l 

OP  0980 

000  2 

l 

OP  96  (DNMI 
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2 

AFIS/INOOC 

6017 

2 
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1 
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1 
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